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The potential role of cosmic topological defects has raised interest in the astrophysical community 
for many years now. In this set of notes, we give an introduction to the subject of cosmic topological 
defects and some of their possible observable signatures. We begin with a review of the basics of 
general defect formation and evolution, we briefly comment on some general features of conducting 
cosmic strings and vorton formation, as well as on the possible role of defects as dark energy, to end 
up with cosmic structure formation from defects and some specific imprints in the cosmic microwave 
background radiation from simulated cosmic strings. A detailed, pedagogical explanation of the 
mechanism underlying the tiny level of polarization discovered in the cosmic microwave background 
by the DASI collaboration (and recently confirmed by WMAP) is also given, and a first rough 
' comparison with some predictions from defects is provided. 
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I. INTRODUCTION 

On a cold day, ice forms quickly on the surface of a pond. But it does not grow as a smooth, featureless covering. In- 
stead, the water begins to freeze in many places independently, and the growing plates of ice join up in random fashion, 
leaving zig-zag boundaries between them. These irregular margins are an example of what physicists call "topological 
defects" - defects because they are places where the crystal structure of the ice is disrupted, and topological because 
an accurate description of them involves ideas of symmetry embodied in topology, the branch of mathematics that 
02 ■ focuses on the study of continuous surfaces. 

Current theories of particle physics likewise predict that a variety of topological defects would almost certainly 
have formed during the early evolution of the universe. Just as water turns to ice (a phase transition) when the 
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' temperature drops, so the interactions between elementary particles run through distinct phases as the typical energy 
of those particles falls with the expansion of the universe. When conditions favor the appearance of a new phase, it 
generally crops up in many places at the same time, and when separate regions of the new phase run into each other, 
topological defects are the result. The detection of such structures in the modern universe would provide precious 
information on events in the earliest instants after the Big Bang. Their absence, on the other hand, would force a 
major revision of current physical theories. 

The aim of this set of Lectures is to introduce the reader to the subject of cosmology from topological defects. 
We begin with a review of the basics of defect formation and evolution, to get a grasp of the overall picture. We 
will see that defects are generically predicted to exist in most interesting models of high energy physics trying to 
describe the early universe. The basic elements of the standard cosmology, with its successes, shortcomings, and new 
developments, are covered elsewhere in this volume. See for example the lecture notes by Rocky Kolb on Astroparticle 
Physics, Ed Copeland's material on String / M-Theory Cosmology, and Jim Bartlett's Observational Cosmology. So 
we will not devote much space to these topics here. Rather, we will focus on some specific subjects. We will first 
briefly comment on conducting cosmic strings and one of their most important predictions for cosmology, namely, 
the existence of equilibrium configurations of string loops, dubbed vortons. We will then pass on to study some key 
signatures that a network of defects would produce on the cosmic microwave background (CMB) radiation, e.g., the 
CMB bispectrum of the temperature anisotropies from a simulated model of cosmic strings. Miscellaneous topics 
also reviewed below are, for example, the way in which these cosmic entities lead to large-scale structure formation 
and some astrophysical footprints left by the various defects, and we will discuss the possibility of isolating their 
effects by astrophysical observations. Also, we include a short, detailed discussion of CMB polarization and some 
brief comparison with the predictions from cosmic defects. 
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Many areas of modern research directly related to cosmic defects are unfortunately not covered in these notes. 
The subject is now so vast -and beyond the possibilities of a single review- that we suggest the reader to consult 
some of the excellent recent literature already available. So, have a look, for example, to the report by Achucarro 
& Vachaspati [2000] for a treatment of semilocal and electroweak strings 1 , and to [Vachaspati, 2001] for a review of 
certain topological defects, like monopoles, domain walls and, again, electroweak strings, virtually not covered here. 
For conducting defects, cosmic strings in particular, see for example [Gangui & Peter, 1998] for a brief overview of many 
different astrophysical and cosmological phenomena, [Gangui, 2001b], from which this review borrows (and updates) 
a great deal of material, for a treatment of conducting cosmic strings and one of their most important predictions 
for cosmology, namely, the existence of equilibrium configurations of string loops, dubbed vortons. Finally, refer to 
the comprehensive colorful lecture notes by Carter [1997] on the dynamics of branes with applications to conducting 
cosmic strings and vortons. If your are in cosmological structure formation, Durrer [2000] presents a good review 
of modern developments on global topological defects and their relation to CMB anisotropics, while Magueijo & 
Brandcnbcrgcr [2000] give a set of imaginative lectures with an update on local string models of large-scale structure 
formation and also baryogenesis with cosmic defects. Finally, Durrer, Kunz & Melchiorri [2002] give a complete 
update of cosmic structure formation with global defects, including detailed analyses of correlators, mixed models, 
and the resulting matter and CMB power spectra. 

The interdisciplinary subject of topological defects in the cosmos and the lab is nicely covered in the proceedings 
of the school held aux Houches on topological defects and non-equilibrium dynamics, edited by Bunkov & Godfrin 
[2000]; the ensemble of lectures in this volume, together with the recent review by Kibble [2002], give an exhaustive 
illustration of this fast developing area of research, which includes various fields of physics, like low-temperature 
condensed-matter, liquid crystals, astrophysics and high-energy physics. Finally, all of the above can also be found 
in the concise review by Hindmarsh & Kibble [1995], particularly concerned with the physics and cosmology of cosmic 
strings, and in the monograph by Vilenkin & Shellard [2000] on cosmic strings and other topological defects. 

A. How defects form 

A central concept of particle physics theories attempting to unify all the fundamental interactions is the concept 
of symmetry breaking. As the universe expanded and cooled, first the gravitational interaction, and subsequently all 
other known forces would have begun adopting their own identities. In the context of the standard hot Big Bang 
theory the spontaneous breaking of fundamental symmetries is realized as a phase transition in the early universe. 
Such phase transitions have several exciting cosmological consequences and thus provide an important link between 
particle physics and cosmology. 

There are several symmetries which are expected to break down in the course of time. In each of these transitions 
the space-time gets 'oriented' by the presence of a hypothetical force field called the 'Higgs field', named for Peter 
Higgs, pervading all the space. This field orientation signals the transition from a state of higher symmetry to a final 
state where the system under consideration obeys a smaller group of symmetry rules. As an every-day analogy we 
may consider the transition from liquid water to ice; the formation of the crystal structure ice (where water molecules 
are arranged in a well defined lattice), breaks the symmetry possessed when the system was in the higher temperature 
liquid phase, when every direction in the system was equivalent. In the same way, it is precisely the orientation in 
the Higgs field which breaks the highly symmetric state between particles and forces. 

Having built a model of elementary particles and forces, particle physicists and cosmologists are today embarked on 
a difficult search for a theory that unifies all the fundamental interactions. As we mentioned, an essential ingredient 
in all major candidate theories is the concept of symmetry breaking. Experiments have determined that there are four 
physical forces in nature; in addition to gravity these are called the strong, weak and electromagnetic forces. Close to 
the singularity of the hot Big Bang, when energies were at their highest, it is believed that these forces were unified 
in a single, all-encompassing interaction. As the universe expanded and cooled, first the gravitational interaction, 
then the strong interaction, and lastly the weak and the electromagnetic forces would have broken out of the unified 
scheme and adopted their present distinct identities in a series of symmetry breakings. 

Theoretical physicists are still struggling to understand how gravity can be united with the other interactions, but for 
the unification of the strong, weak and electromagnetic forces plausible theories exist. Indeed, force-carrying particles 
whose existence demonstrated the fundamental unification of the weak and electromagnetic forces into a primordial 
"electroweak" force - the W and Z bosons - were discovered at CERN, the European accelerator laboratory, in 1983. 
In the context of the standard Big Bang theory, cosmological phase transitions are produced by the spontaneous 



1 Animations of semilocal and electroweak string formation and evolution can be found at http://www.nersc.gov/-borrill/ 
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breaking of a fundamental symmetry, such as the electroweak force, as the universe cools. For example, the electroweak 
interaction broke into the separate weak and electromagnetic forces when the observable universe was 10 -12 seconds 
old, had a temperature of 10 15 degrees Kelvin, and was only one part in 10 15 of its present size. There are also other 
phase transitions besides those associated with the emergence of the distinct forces. The quark-hadron confinement 
transition, for example, took place when the universe was about a microsecond old. Before this transition, quarks 
- the particles that would become the constituents of the atomic nucleus - moved as free particles; afterward, they 
became forever bound up in protons, neutrons, mesons and other composite particles. 

As we said, the standard mechanism for breaking a symmetry involves the hypothetical Higgs field that pervades 
all space. As the universe cools, the Higgs field can adopt different ground states, also referred to as different vacuum 
states of the theory. In a symmetric ground state, the Higgs field is zero everywhere. Symmetry breaks when the 
Higgs field takes on a finite value (see Figure QJ. 
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FIG. 1: Temperature-dependent effective potential for a first-order phase transition for the Higgs held. For very high 
temperatures, well above the critical one T c , the potential possesses just one minimum for the vanishing value of the Higgs 
held. Then, when the temperature decreases, a whole set of minima develops (it may be two or more, discrete or continuous, 
depending of the type of symmetry under consideration). Below T c , the value <j> — stops being the global minimum and the 
system will spontaneously choose a new (lower) one, say 4> = rjexp(i9) (for complex <p) for some angle 9 and nonvanishing r\, 
amongst the available ones. This choice signals the breakdown of the symmetry in a cosmic phase transition and the generation 
of random regions of conflicting held orientations 9. In a cosmological setting, the merging of these domains gives rise to cosmic 
defects. 

Kibble [1976] first saw the possibility of defect formation when he realized that in a cooling universe phase transitions 
proceed by the formation of uncorrelated domains that subsequently coalesce, leaving behind relics in the form of 
defects. In the expanding universe, widely separated regions in space have not had enough time to 'communicate' 
amongst themselves and are therefore not correlated, due to a lack of causal contact. It is therefore natural to suppose 
that different regions ended up having arbitrary orientations of the Higgs field and that, when they merged together, 
it was hard for domains with very different preferred directions to adjust themselves and fit smoothly. In the interfaces 
of these domains, defects form. Such relic 'flaws' are unique examples of incredible amounts of energy and this feature 
attracted the minds of many cosmologists. 

B. Phase transitions and finite temperature field theory 

Phase transitions are known to occur in the early universe. Examples we mentioned are the quark to hadron 
(confinement) transition, which QCD predicts at an energy around I GeV, and the electroweak phase transition at 
about 250 GeV. Within grand unified theories (GUT), aiming to describe the physics beyond the standard model, 
other phase transitions are predicted to occur at energies of order 10 15 GeV; during these, the Higgs field tends to 
fall towards the minima of its potential while the overall temperature of the universe decreases as a consequence of 
the expansion. 

A familiar theory to make a bit more quantitative the above considerations is the A|</>| 4 theory, 

C- = \\d^? + \ml\<t>?-^\<l>W (i) 

with rriQ > 0. The second and third terms on the right hand side yield the usual 'Mexican hat' potential for the 
complex scalar field. For energies much larger than the critical temperature, T c , the fields are in the so-called 'false' 
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vacuum: a highly symmetric state characterized by a vacuum expectation value = 0. But when energies decrease 
the symmetry is spontaneously broken: a new 'true' vacuum develops and the scalar field rolls down the potential and 

sits onto one of the degenerate new minima. In this situation the vacuum expectation value becomes ([</>[) = 6m 2 ,/ A. 

Research done in the 1970's in finite-temperature field theory [Weinberg, 1974; Dolan & Jackiw, 1974; Kirzhnits & 
Linde, 1974] has led to the result that the temperature-dependent effective potential can be written down as 

W(|0|) = -im 2 (T)|^| 2 + A|^ (2) 

with T 2 = 24m§/A, m 2 (T) = m§(l - T 2 /T 2 ), and (\(j)\) 2 = 6m 2 (T)/A. We easily see that when T approaches T c from 
below the symmetry is restored, and again we have (\4>\) = 0. In condensed-matter jargon, the transition described 
above is second-order [Mermin, 1979]. 2 

C. The Kibble mechanism 

The model described in the last subsection is an example in which the transition may be second-order. As we saw, 
for temperatures much larger than the critical one the vacuum expectation value of the scalar field vanishes at all 
points of space, whereas for T < T c it evolves smoothly in time towards a non vanishing Both thermal and 

quantum fluctuations influence the new value taken by an d therefore it has no reasons to be uniform in space. 
This leads to the existence of domains wherein the (|</>(x)|) is coherent and regions where it is not. The consequences 
of this fact are the subject of this subsection. 

Phase transitions can also be first-order proceeding via bubble nucleation. At very high energies the symmetry 
breaking potential has (\4>\) = as the only vacuum state. When the temperature goes down to T c a set of vacua, 
degenerate to the previous one, develops. However this time the transition is not smooth as before, for a potential 
barrier separates the old (false) and the new (true) vacua (see, e.g. Figure Provided the barrier at this small 
temperature is high enough, compared to the thermal energy present in the system, the field 4> will remain trapped 
in the false vacuum state even for small (< T c ) temperatures. Classically, this is the complete picture. However, 
quantum tunneling effects can liberate the field from the old vacuum state, at least in some regions of space: there is 
a probability per unit time and volume in space that at a point x a bubble of true vacuum will nucleate. The result is 
thus the formation of bubbles of true vacuum with the value of the field in each bubble being independent of the value 
of the field in all other bubbles. This leads again to the formation of domains where the fields are correlated, whereas 
no correlation exits between fields belonging to different domains. Then, after creation the bubble will expand at the 
speed of light surrounded by a 'sea' of false vacuum domains. As opposed to second-order phase transitions, here the 
nucleation process is extremely inhomogeneous and (|0(:e)|) is not a continuous function of time. 

Let us turn now to the study of correlation lengths and their role in the formation of topological defects. One 
important feature in determining the size of the domains where (|0(a;)|) is coherent is given by the spatial correlation 
of the field cf). Simple field theoretic considerations [see, e.g., Copeland, 1993] for long wavelength fluctuations of <f> 
lead to different functional behaviors for the correlation function G(r) = ^</>(ri)0(r2)), where we noted r = \r% — r%\. 
What is found depends radically on whether the wanted correlation is computed between points in space separated 
by a distance r much smaller or much larger than a characteristic length = m(T) ~ VA |(^)|, known as the 

correlation length. Then, we have G(r) ~ ^ exp(— |) for r » £ , while G(r) ~ ^ for r << £ . 

This tells us that domains of size £ ~ m _1 arise where the field <f> 1S correlated. On the other hand, well beyond £ 
no correlations exist and thus points separated apart by r >> £ will belong to domains with in principle arbitrarily 
different orientations of the Higgs field. This in turn leads, after the merging of these domains in a cosmological 
setting, to the existence of defects, where field configurations fail to match smoothly. 

However, when T — > T c we have m — > and so £ — > oo, suggesting perhaps that for all points of space the field 
4> becomes correlated. This fact clearly violates causality. The existence of particle horizons in cosmological models 
(proportional to the inverse of the Hubble parameter H _1 ) constrains microphysical interactions over distances beyond 
this causal domain. Therefore we get an upper bound to the correlation length as £ < H^ 1 ~ t. 

The general feature of the existence of uncorrelated domains has become known as the Kibble mechanism [Kibble, 
1976] and it seems to be generic to most types of phase transitions. 



2 In a first-order phase transition the order parameter (e.g., < \(j>\ > in our case) is not continuous. It may proceed by bubble nucleation 
[Callan & Coleman, 1977; Linde, 1983b] or by spinoidal decomposition [Langer, 1992]. Phase transitions can also be continuous 
second-order processes. The 'order' depends sensitively on the ratio of the coupling constants appearing in the Lagrangian. 
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D. A survey of topological defects 

Different models for the Higgs field lead to the formation of a whole variety of topological defects, with very 
different characteristics and dimensions. Some of the proposed theories have symmetry breaking patterns leading 
to the formation of 'domain walls' (mirror reflection discrete symmetry): incredibly thin planar surfaces trapping 
enormous concentrations of mass-energy which separate domains of conflicting field orientations, similar to two- 
dimensional sheet - like structures found in fcrromagncts. Within other theories, cosmological fields get distributed in 
such a way that the old (symmetric) phase gets confined into a finite region of space surrounded completely by the 
new (non-symmetric) phase. This situation leads to the generation of defects with linear geometry called 'cosmic 
strings'. Theoretical reasons suggest these strings (vortex lines) do not have any loose ends in order that the two 
phases not get mixed up. This leaves infinite strings and closed loops as the only possible alternatives for these defects 
to manifest themselves in the early universe 3 . 




FIG. 2: In a simple model of symmetry breaking, the initial symmetric ground state of the Higgs field (yellow central dot) 
can fall into the left- or right-hand valley of a double-well energy potential (light and dark dots). In a cosmic phase transition, 
regions of the new phase appear randomly and begin to grow and eventually merge as the transition proceeds toward completion 
(middle). Regions in which the symmetry has broken the same way can coalesce, but where regions that have made opposite 
choices encounter each other, a topological defect known as a domain wall forms (right). Across the wall, the Higgs field has to 
go from one of the valleys to the other (in the left panel), and must therefore traverse the energy peak. This creates a narrow 
planar region of very high energy, in which the symmetry is locally unbroken. 

With a bit more abstraction scientists have even conceived other (semi) topological defects, called 'textures'. These 
are conceptually simple objects, yet, it is not so easy to imagine them for they are just global field configurations living 
on a three-sphere vacuum manifold (the minima of the effective potential energy), whose non linear evolution perturbs 
spacetime. Turok [1989] was the first to realize that many unified theories predicted the existence of peculiar Higgs 
field configurations known as (texture) knots, and that these could be of potential interest for cosmology. Several 
features make these defects interesting. In contrast to domain walls and cosmic strings, textures have no core and 
thus the energy is more evenly distributed over space. Secondly, they are unstable to collapse and it is precisely 
this last feature which makes these objects cosmologically relevant, for this instability makes texture knots shrink to 
a microscopic size, unwind and radiate away all their energy. In so doing, they generate a gravitational field that 
perturbs the surrounding matter in a way which can seed structure formation. 

E. Conditions for their existence: topological criteria 

Let us now explore the conditions for the existence of topological defects. It is widely accepted that the final 
goal of particle physics is to provide a unified gauge theory comprising strong, weak and electromagnetic interactions 
(and some day also gravitation). This unified theory is to describe the physics at very high temperatures, when 
the age of the universe was slightly bigger than the Planck time. At this stage, the universe was in a state with the 



'Monopole' is another possible topological defect; we defer its discussion to the next subsection. Cosmic strings bounded by monopoles 
is yet another possibility in GUT phase transitions of the kind, e.g., G — » K X t/(l) — » K. The first transition yields monopoles carrying 
a magnetic charge of the U(l) gauge field, while in the second transition the magnetic field in squeezed into flux tubes connecting 
monopoles and antimonopoles [Langacker & Pi, 1980]. 
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-ko(M) M disconnected Domain Walls 

■ki(M) non contractible loops in M Cosmic Strings 

■K2{M) non contractible 2-spheres in M. Monopoles 

■k?,{M) non contractible 3-spheres in M Textures 



TABLE I: The topology of M determines the type of defect that will arise. 

highest possible symmetry, described by a symmetry group G, and the Lagrangian modeling the system of all possible 
particles and interactions present should be invariant under the action of the elements of G. 

As we explained before, the form of the finite temperature effective potential of the system is subject to variations 
during the cooling down evolution of the universe. This leads to a chain of phase transitions whereby some of the 
symmetries present in the beginning are not present anymore at lower temperatures. The first of these transitions 
may be described as G— >H, where now H stands for the new (smaller) unbroken symmetry group ruling the system. 
This chain of symmetry breakdowns eventually ends up with SU(3)xSU(2)xU(l), the symmetry group underlying 
the 'standard model' of particle physics. 

A broken symmetry system (with a Mexican-hat potential for the Higgs field) may have many different minima 
(with the same energy), all related by the underlying symmetry. Passing from one minimum to another is included 
as one of the symmetries of the original group G, and the system will not change due to one such transformation. 
If a certain field configuration yields the lowest energy state of the system, transformations of this configuration by 
the elements of the symmetry group will also give the lowest energy state. For example, if a spherically symmetric 
system has a certain lowest energy value, this value will not change if the system is rotated. 

The system will try to minimize its energy and will spontaneously choose one amongst the available minima. Once 
this is done and the phase transition achieved, the system is no longer ruled by G but by the symmetries of the smaller 
group H. So, if G^H and the system is in one of the lowest energy states (call it Si), transformations of S\ to 5*2 
by elements of G will leave the energy unchanged. However, transformations of S\ by elements of H will leave S\ 
itself (and not just the energy) unchanged. The many distinct ground states of the system Si, 62, • • • are given by all 
transformations of G that are not related by elements in H. This space of distinct ground states is called the vacuum 
manifold and denoted M. . So, M. is the space of all elements of G in which elements related by transformations in H 
have been identified. Mathematicians call it the coset space and denote it G/H. We then have M. = G/H. 

The importance of the study of the vacuum manifold lies in the fact that it is precisely the topology of M. what 
determines the type of defect that will arise. Homotopy theory tells us how to map M. into physical space in a 
non-trivial way, and what ensuing defect will be produced. For instance, the existence of non contractible loops in 
M. is the requisite for the formation of cosmic strings. In formal language this comes about whenever we have the 
first homotopy group iti(M) ^ 1, where 1 corresponds to the trivial group. If the vacuum manifold is disconnected 
we then have ttq(M) ^ 1, and domain walls are predicted to form in the boundary of these regions where the field <p 
is away from the minimum of the potential. Analogously, if ^(.M) 7^ 1 it follows that the vacuum manifold contains 
non contractible two-spheres, and the ensuing defect is a monopole. Textures arise when M. contains non contractible 
three-spheres and in this case it is the third homotopy group, 713 (M), the one that is non trivial. We summarize this 
in Table [Q. 

II. DEFECTS IN THE UNIVERSE 

Generically topological defects will be produced if the conditions for their existence are met. Then for example if 
the unbroken group H contains a disconnected part, like an explicit U(l) factor (something that is quite common in 
many phase transition schemes discussed in the literature), monopoles will be left as relics of the transition. This is 
due to the fundamental theorem on the second homotopy group of coset spaces [Mermin, 1979], which states that for 
a simply-connected covering group G we have 4 

7r 2 (G/H)S*7n(Ho), (3) 



4 The isomorfism between two groups is noted as =. Note that by using the theorem we therefore can reduce the computation of n2 for a 
coset space to the computation of tti for a group. A word of warning: the focus here is on the physics and the mathematically-oriented 
reader should bear this in mind, especially when we will become a bit sloppy with the notation. In case this happens, consult the book 
[Steenrod, 1951] for a clear exposition of these matters. 
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with Ho being the component of the unbroken group connected to the identity. Then we see that since monopolcs 
are associated with unshrinkable surfaces in G/H, the previous equation implies their existence if H is multiply- 
connected. The reader may guess what the consequences are for GUT phase transitions: in grand unified theories a 
semi-simple gauge group G is broken in several stages down to H = SU(3) xU(l). Since in this case 7Ti(H) = Z, the 
integers, we have 7r 2 (G/H) ^ 1 and therefore gauge monopole solutions exist [Preskill, 1979]. 

A. Local and global monopoles and domain walls 

Monopoles are yet another example of stable topological defects. Their formation stems from the fact that the 
vacuum expectation value of the symmetry breaking Higgs field has random orientations (((j) a ) pointing in different 
directions in group space) on scales greater than the horizon. One expects therefore to have a probability of order 
unity that a monopole configuration will result after the phase transition (cf. the Kibble mechanism). Thus, about 
one monopole per Hubble volume should arise and we have for the number density n monop <~ l/H~ 3 <~ T^/mp, where 
T c is the critical temperature and nip is Planck mass, when the transition occurs. We also know the entropy density 
at this temperature, s ~ T 3 , and so the monopole to entropy ratio is n monop /s ~ 100(T c /mp) 3 . In the absence 
of non-adiabatic processes after monopole creation this constant ratio determines their present abundance. For the 
typical value T c ~ 10 14 GeV we have n monop /s ~ 1CP 13 . This estimate leads to a present ^l m ono P h 2 — 10 11 , for the 
superheavy monopoles m monop ~ 10 16 GeV that are created 5 . This value contradicts standard cosmology and the 
presently most attractive way out seems to be to allow for an early period of inflation: the massive entropy production 
will hence lead to an exponential decrease of the initial n monop /s ratio, yielding £l m£mop consistent with observations. 6 
In summary, the broad-brush picture one has in mind is that of a mechanism that could solve the monopole problem 
by 'weeping' these unwanted relics out of our sight, to scales much bigger than the one that will eventually become 
our present horizon today. 

Note that these arguments do not apply for global monopoles as these (in the absence of gauge fields) possess 
long-range forces that lead to a decrease of their number in comoving coordinates. The large attractive force between 
global monopoles and antimonopoles leads to a high annihilation probability and hence monopole over-production 
does not take place. Simulations performed by Bennett & Rhie [1990] showed that global monopole evolution rapidly 
settles into a scale invariant regime with only a few monopoles per horizon volume at all times. 

Given that global monopoles do not represent a danger for cosmology one may proceed in studying their observable 
consequences. The gravitational fields of global monopoles may lead to matter clustering and CMB anisotropics. 
Given an average number of monopoles per horizon of <~ 4, Bennett & Rhie [1990] estimate a scale invariant spectrum 
of fluctuations (dp/ ' p)n ~ 30G?] 2 at horizon crossing 7 . In a subsequent paper they simulate the large-scale CMB 
anisotropies and, upon normalization with COBE-DMR, they get roughly Grj 2 ~6x 10~ 7 in agreement with a GUT 
energy scale n [Bennett & Rhie, 1993]. However, as we will see in the CMB sections below, current estimates for the 
angular power spectrum of global defects do not match the most recent observations, their main problem being the 
lack of power on the degree angular scale once the spectrum is normalized to COBE on large scales [Durrer et al., 
1996; Durrer et al., 2002]. 

Let us concentrate now on domain walls, and briefly try to show why they are not welcome in any cosmological 
context [at least in the simple version we here consider - there is always room for more complicated (and contrived) 
models]. If the symmetry breaking pattern is appropriate at least one domain wall per horizon volume will be formed. 
The mass per unit surface of these two-dimensional objects is given by ~ A 1 / 2 /; 3 , where A as usual is the coupling 
constant in the symmetry breaking potential for the Higgs field. Domain walls are generally horizon-sized and therefore 
their mass is given by ~ \ 1 / 2 if' H~ 2 . This implies a mass energy density roughly given by pdw ~ rpt^ 1 and we may 
readily see now how the problem arises: the critical density goes as p cr it ~ t~ 2 which implies Q.Dw(t) ~ (i)/mp) 2 i)(. 
Taking a typical GUT value for r\ we get Clowit ~ 10~ 35 sec) ~ 1 already at the time of the phase transition. It is 
not hard to imagine that today this will be at variance with observations; in fact we get flowit ~ 10 18 sec) ~ 10 52 . 



5 These are the actual figures for a gauge SU(5) GUT second-order phase transition. Preskill [1979] has shown that in this case monopole 
antimonopolc annihilation is not effective to reduce their abundance. Guth & Weinberg [1983] did the case for a first-order phase 
transition and drew qualitatively similar conclusions regarding the excess of monopoles. 

6 The inflationary expansion reaches an end in the so-called reheating process, when the enormous vacuum energy driving inflation is 
transferred to coherent oscillations of the inflaton field. These oscillations will in turn be damped by the creation of light particles (e.g., 
via preheating) whose final fate is to thermalise and reheat the universe. 

7 The spectrum of density fluctuations on smaller scales has also been computed. They normalize the spectrum at 8h~ 1 Mpc and 
agreement with observations lead them to assume that galaxies are clustered more strongly than the overall mass density, this implying 
a 'biasing' of a few [see Bennett, Rhie & Weinberg, 1993 for details]. 
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This indicates that models where domain walls are produced are tightly constrained, and the general feeling is that 
it is best to avoid them altogether [see Kolb & Turner, 1990 for further details; see also Dvali et al., 1998, Pogosian 
& Vachaspati, 2000 8 and Alexander et al., 1999 for an alternative solution]. 



B. Are defects inflated away? 



It is important to realize the relevance that the Kibble's mechanism has for cosmology; nearly every sensible grand 
unified theory (with its own symmetry breaking pattern) predicts the existence of defects. We know that an early era 
of inflation helps in getting rid of the unwanted relics. One could well wonder if the very same Higgs field responsible 
for breaking the symmetry would not be the same one responsible for driving an era of inflation, thereby diluting 
the density of the relic defects. This would get rid not only of (the unwanted) monopoles and domain walls but also 
of any other (cosmologically appealing) defect. Let us follow [Brandenbcrger, 1993] and sketch why this actually 
does not occur. Take first the symmetry breaking potential of Eq. Q at zero temperature and add to it a harmless 
^-independent term 3m /(2A). This will not affect the dynamics at all. Then we are led to 

VW = ± (<t> 2 -r, 2 ) 2 , (4) 

with r\ = (<om 2 / X) 1 / 2 the symmetry breaking energy scale, and where for the present heuristic digression we just took 
a real Higgs field. Consider now the equation of motion for </>, 

dV A , q , , , , 

= --0 3 +m 2 0«m 2 0, (5) 

for <p « rj very near the false vacuum of the effective Mexican hat potential and where, for simplicity, the expansion 
of the universe and possible interactions of cf> with other fields were neglected. The typical time scale of the solution 
is t ~ i7i~ . For an inflationary epoch to be effective we need r >> H -1 , i.e., a sufficiently large number of e-folds 
of slow-rolling solution. Note, however, that after some e-folds of exponential expansion the curvature term in the 
Friedmann equation becomes subdominant and we have H 2 ~ 8irG V(0)/3 ~ (2irm 2 /3)(i]/mp) 2 . So, unless rj > mp, 
which seems unlikely for a GUT phase transition, we are led to r << H^ 1 and therefore the amount of inflation is not 
enough for getting rid of the defects generated during the transition by hiding them well beyond our present horizon. 

Recently, there has been a large amount of work in getting defects, particularly cosmic strings, after post-inflationary 
preheating. Reaching the latest stages of the inflationary phase, the inflaton field oscillates about the minimum of 
its potential. In doing so, parametric resonance may transfer a huge amount of energy to other fields leading to 
cosmologically interesting nonthermal phase transitions. Just like thermal fluctuations can restore broken symmetries, 
here also, these large fluctuations may lead to the whole process of defect formation again. Numerical simulations 
employing potentials similar to that of Eq. have shown that strings indeed arise for values 77 ~ 10 16 GeV [Tkachev 
et al., 1998, Kasuya & Kawasaki, 1998]. Hence, preheating after inflation helps in generating cosmic defects. 



C. Cosmic strings 



Cosmic strings are without any doubt the topological defect most thoroughly studied, both in cosmology and 
solid-state physics (vortices). The canonical example, also describing flux tubes in superconductors, is given by the 
Lagrangian 




with — dtpA,,], where A v is the gauge field and the covariant derivative is = + ieA^, with e the gauge 
coupling constant. This Lagrangian is invariant under the action of the Abelian group G = U(l), and the spontaneous 
breakdown of the symmetry leads to a vacuum manifold Ai that is a circle, S , i.e., the potential is minimized for 
4> = r]exp(i6), with arbitrary < 9 < 2ir. Each possible value of 9 corresponds to a particular 'direction' in the field 
space. 



Animations of monopoles colliding with domain walls can be found in 'LEP' page at http://theory.ic.ac.uk/-LEP/figures.html 
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FIG. 3: The complex scalar Higgs field evolves in a temperature-dependent potential V(<j>). At high temperatures (violet narrow 
surface) the vacuum expectation value of the field lies at the bottom of V . For lower temperatures, the potential adopts the 
"Mexican hat" form (yellow surface) and the held spontaneously chooses one amongst the new available (degenerate) lowest 
energy states (the violet circle along the valley of the hat). This isolates a single value/direction for the phase of the Held, 
spontaneously breaking the symmetry possessed by the system at high energies. Different regions of the universe, with no 
causal connection, will end up having arbitrarily different directions for the held (arrows on the right). As separate regions 
of broken symmetry merge, it is not always possible for the field orientations to match. It may happen that a closed loop in 
physical space intersects regions where the Higgs phase varies from to 2n (red arrows, corresponding to the red dashed-line 
on the left panel). In that situation, a cosmic string will pass somewhere inside the loop. On the contrary, green arrows (and 
green dashed-line on the left panel) show a situation where no string is formed after the phase transition. 

Now, as we have seen earlier, due to the overall cooling down of the universe, there will be regions where the scalar 
field rolls down to different vacuum states. The choice of the vacuum is totally independent for regions separated 
apart by one correlation length or more, thus leading to the formation of domains of size £ ~ V • When these 
domains coalesce they give rise to edges in the interface. If we now draw a imaginary circle around one of these edges 
and the angle 9 varies by 2ir then by contracting this loop we reach a point where we cannot go any further without 
leaving the manifold M. This is a small region where the variable 9 is not defined and, by continuity, the field should 
be <p — 0. In order to minimize the spatial gradient energy these small regions line up and form a line-like defect 
called cosmic string. 

The width of the string is roughly m^ 1 ~ (y/Xr)) -1 , m$ being the Higgs mass. The string mass per unit length, 
or tension, is /i ~ rj 2 . This means that for GUT cosmic strings, where rj ~ 10 16 GeV, we have G/i ~ 10 -6 . We will 
see below that the dimensionless combination G/i, present in all signatures due to strings, is of the right order of 
magnitude for rendering these defects cosmologically interesting. 

There is an important difference between global and gauge (or local) cosmic strings: local strings have their energy 
confined mainly in a thin core, due to the presence of gauge fields that cancel the gradients of the field outside 
of it. Also these gauge fields make it possible for the string to have a quantized magnetic flux along the core. On 
the other hand, if the string was generated from the breakdown of a global symmetry there are no gauge fields, just 
Goldstone bosons, which, being massless, give rise to long-range forces. No gauge fields can compensate the gradients 
of 4> this time and therefore there is an infinite string mass per unit length. 

Just to get a rough idea of the kind of models studied in the literature, consider the case G = 50(10) that is broken 
to H = SU(5) x 2>%, For this pattern we have 7Ti(A4) = Z2, which is clearly non trivial and therefore cosmic strings 
are formed [Kibble et al, 1982]. 9 

D. String loops and scaling 

We saw before the reasons why gauge monopoles and domain walls were a bit of a problem for cosmology. Essentially, 
the problem was that their energy density decreases more slowly than the critical density with the expansion of the 



In the analysis one uses the fundamental theorem stating that, for a simply— connected Lie group G breaking down to H, we have 
7Ti(G/H) Stto(H); see [Hilton, 1953]. 
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universe. This fact resulted in their contribution to fidcf (the density in defects normalized by the critical density) 
being largely in excess compared to 1, hence in blatant conflict with modern observations. The question now arises as 
to whether the same might happened with cosmic strings. Are strings dominating the energy density of the universe? 
Fortunately, the answer to this question is no; strings evolve in such a way to make their density ^strings °c rj 2 t~ 2 . Hence, 
one gets the same temporal behavior as for the critical density. The result is that firings ~ Gfi ~ (rj/mp) 2 ~ 10~ 6 
for GUT strings, i.e., we get an interestingly small enough, constant fraction of the critical density of the universe 
and strings never upset standard observational cosmology. 



3 - DIMENSIONAL 

PHYSICAL SPACE PHYSICAL SPACE 




FIG. 4: We can now extend the mechanism shown in the previous figure to the full three-dimensional space. Regions of the 
various planes that were traversed by strings can be superposed to show the actual location of the cosmic string (left panel). 
The figure on the right panel shows why we are sure a string crosses the plane inside the loop in physical space (the case with 
red arrows in the previous figure). Continuity of the field imposes that if we gradually contract this loop the direction of the 
field will be forced to wind "faster". In the limit in which the loop reduces to a point, the phase is no longer deBned and the 
vacuum expectation value of the Higgs held has to vanish. This corresponds to the central tip of the Mexican hat potential 
in the previous figure and is precisely the locus of the false vacuum. Cosmic strings are just that, narrow, extremely massive 
line-like regions in physical space where the Higgs held adopts its high-energy false vacuum state. 



Now, why this is so? The answer is simply the efficient way in which a network of strings looses energy. The 
evolution of the string network is highly nontrivial and loops are continuously chopped off from the main infinite 
strings as the result of (self) intersections within the infinite-string network. Once they are produced, loops oscillate 
due to their huge tension and slowly decay by emitting gravitational radiation. Thus, energy is transferred from the 
cosmic string network to radiation. 10 

It turns out from simulations that most of the energy in the string network (roughly a 80%) is in the form of infinite 
strings. Soon after formation one would expect long strings to have the form of random- walk with characteristic step 
given by the correlation length £. Also, the typical distance between long string segments should also be of order 
£. Monte Carlo simulations show that these strings are Brownian on sufficiently large scales, which means that the 
length 1 of a string is related to the end-to-end distance d of two given points along the string (with d>$) in the 
form 



i = d 2 /e (7) 

What remains of the energy is given in the form of closed loops with no preferred length scale (a scale invariant 
distribution) which implies that the number density of loops having sizes between R and R + dR follows just from 
dimensional analysis 

an-ioops oc --j (8) 

which is just another way of saying that ni oops cx 1/R 3 , loops behave like normal nonrelativistic matter. The actual 
coefficient, as usual, comes from string simulations. 

There are both analytical and numerical indications in favor of the existence of a stable "scaling solution" for the 
cosmic string network. After generation, the network quickly evolves in a self similar manner with just a few infinite 



High-resolution cosmic string simulations can be found in the Cambridge cosmology page at 
http : //www. damtp . cam. ac .uk/user/gr/public/cs_evol .html 
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string segments per Hubble volume and Hubble time. A heuristic argument for the scaling solution due to Vilenkin 
[1985] is as follows. 

If we take v(t) to be the mean number of infinite string segments per Hubble volume, then the energy density in 
infinite strings brings = Ps is 

Ps {t) = v{i)^r 2 = V {t)^r 2 . (9) 

Now, v strings will typically have v intersections, and so the number of loops ni oops (t) = ni(t) produced per unit 
volume will be proportional to v 2 . We find 

dm ~ v 2 TC^dR. (10) 
Hence, recalling now that the loop sizes grow with the expansion like R oc t we have 

dn\(t) 



dt 



pv 2 r 4 (11) 



where p is the probability of loop formation per intersection, a quantity related to the intercommuting probability, 
both roughly of order 1. We are now in a position to write an energy conservation equation for strings plus loops in 
the expanding universe. Here it is 

dp s , 3 dn\ dn\ 

+ » "A* (12) 

where m\ — fit is just the loop mass and where the second on the left hand side is the dilution term 3Hp s for an 
expanding radiation-dominated universe. The term on the right hand side amounts to the loss of energy from the 
long string network by the generation of small closed loops. Plugging Eqs. © and (|llfl into (|12|) Vilenkin finds the 
following kinetic equation for v{t) 

dv v v 2 

ti-2t~- p T (13) 

with p ~ 1. Thus if v 3> 1 then dvjdt < and v tends to decrease in time, while if v <C 1 then dvjdt > and v 
increases. Hence, there will be a stable solution with v ~ a few. 



E. Global textures 



Whenever a global non-Abelian symmetry is spontaneously and completely broken (e.g. at a grand unification 
scale), global defects called textures are generated. Theories where this global symmetry is only partially broken do 
not lead to global textures, but instead to global monopoles and non-topological textures. As we already mentioned 
global monopoles do not suffer the same constraints as their gauge counterparts: essentially, having no associated 
gauge fields, the long-range forces between pairs of monopoles lead to the annihilation of their eventual excess and 
as a result monopoles scale with the expansion. On the other hand, non-topological textures are a generalization 
that allows the broken subgroup H to contain non-Abelian factors. It is then possible to have tt^ trivial as in, e.g., 
SO(5)— >SO(4) broken by a vector, for which case we have M. — S , the four-sphere [Turok, 1989]. Having explained 
this, let us concentrate in global topological textures from now on. 

Textures, unlike monopoles or cosmic strings, are not well localized in space. This is due to the fact that the 
field remains in the vacuum everywhere, in contrast to what happens for other defects, where the field leaves the 
vacuum manifold precisely where the defect core is. Since textures do not possess a core, all the energy of the field 
configuration is in the form of field gradients. This fact is what makes them interesting objects only when coming 
from global theories: the presence of gauge fields could (by a suitable reorientation) compensate the gradients of 
(f> and yield D^(f> = 0, hence canceling out (gauging away) the energy of the configuration 11 . 



11 This does not imply, however, that the classical dynamics of a gauge texture is trivial. The evolution of the <j>— system will be 
determined by the competing tendencies of the global field to unwind and of the gauge field to compensate the <j> gradients. The result 
depends on the characteristic size L of the texture: in the range m^ 1 << L << m^ 1 ~ (er?) -1 the behavior of the gauge texture 
resembles that of the global texture, as it should, since in the limit very small (e — » 0) the gauge texture turns into a global one 
[Turok & Zadrozny, 1990]. 
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FIG. 5: Global string interactions leading to loop formation. Whenever two string segments intersect, they reconnect or 
intercommute (green and red strings - upper part of the figure). Analogously, if a string intersects itself, it can break off 
a closed loop (green string - bottom part of the figure). In both cases, the interacting string segments Erst suffer a slight 
deformation (due to the long-range forces present for global strings), they subsequently fuse and finally exchange partners. A 
ephemeral unstable amount of energy in the form of a small loop remains in the middle where the energy is high enough to 
place the Higgs field in the false vacuum. It then quickly collapses, radiating away its energy. The situation is roughly the 
same for local strings, as simulations have shown. 



One feature endowed by textures that really makes these defects peculiar is their being unstable to collapse. The 
initial field configuration is set at the phase transition, when <p develops a nonzero vacuum expectation value. <f> lives 
in the vacuum manifold M and winds around M in a non-trivial way on scales greater than the correlation length, 
£ t. The evolution is determined by the nonlinear dynamics of <\>. When the typical size of the defect becomes of the 
order of the horizon, it collapses on itself. The collapse continues until eventually the size of the defect becomes of the 
order of r\~ x , and at that point the energy in gradients is large enough to raise the field from its vacuum state. This 
makes the defect unwind, leaving behind a trivial field configuration. As a result £ grows to about the horizon scale, 
and then keeps growing with it. As still larger scales come across the horizon, knots are constantly formed, since the 
field <f> points in different directions on M. in different Hubble volumes. This is the scaling regime for textures, and 
when it holds simulations show that one should expect to find of order 0.04 unwinding collapses per horizon volume 
per Hubble time [Turok, 1989]. However, unwinding events are not the most frequent feature [Borrill et al, 1994], and 
when one considers random field configurations without an unwinding event the number raises to about 1 collapse 
per horizon volume per Hubble time. 



F. Evolution of global textures 



We mentioned earlier that the breakdown of any non-Abelian global symmetry led to the formation of textures. 
The simplest possible example involves the breakdown of a global SU(2) by a complex doublet (f> a , where the latter 
may be expressed as a four-component scalar field, i.e., a = 1 ... 4. We may write the Lagrangian of the theory much 
in the same way as it was done in Eq. JfjJ, but now we drop the gauge fields (thus the covariant derivatives become 

partial derivatives). Let us take the symmetry breaking potential as follows, V{<j>) — j (|0| 2 — r/ 2 ) . The situation in 
which a global SU(2) in broken by a complex doublet with this potential V is equivalent to the theory where SO(4) is 
broken by a four-component vector to SO (3), by making 4> a take on a vacuum expectation value. We then have the 
vacuum manifold M given by SO(4)/SO(3) = S 3 , namely, a three-sphere with </> a Q = rj 2 . As ^(S 3 ) ^ 1 (in fact, 
7T3(S' 3 ) = Z) we see we will have non-trivial solutions of the field 4> a and global textures will arise. 
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As usual, variation of the action with respect to the field <fi a yields the equation of motion 

^ + 2 fy_vy = _a a X ) (14) 
a oq>° 

where primes denote derivatives with respect to conformal time and V is computed in comoving coordinates. When 
the symmetry in broken three of the initially four degrees of freedom go into massless Goldstone bosons associated 
with the three directions tangential to the vacuum three-sphere. The 'radial' massive mode that remains (m^ ~ sTXf]) 
will not be excited, provided we concentrate on length scales much larger than ml . 

To solve for the dynamics of the field (f> b , two different approaches have been implemented in the literature. The 
first one faces directly the full equation i|14|) . trying to solve it numerically. The alternative to this exploits the fact 
that, at temperatures smaller than T c , the field is constrained to live in the true vacuum. By implementing this fact 
via a Lagrange multiplier 12 we get 

VV = -^^ ; ^ = v\ (is) 

with V M the covariant derivative operator. Eq. (|15f) represents a non-linear sigma model for the interaction of the 
three massless modes [Rajaraman, 1982]. This last approach is only valid when probing length scales larger than the 
inverse of the mass m^ 1 . As we mentioned before, when this condition is not met the gradients of the field are strong 
enough to make it leave the vacuum manifold and unwind. 

The approach (cf. Eqs. I|15|)) is suitable for analytic inspection. In fact, an exact flat space solution was found 
assuming a spherically symmetric ansatz. This solution represents the collapse and subsequent conversion of a texture 
knot into massless Goldstone bosons, and is known as the spherically symmetric self-similar (SSSS) exact unwinding 
solution. We will say no more here with regard to the this solution, but just refer the interested reader to the original 
articles [see, e.g., Turok & Spergel, 1990; Notzold, 1991]. Simulations taking full account of the energy stored in 
gradients of the field, and not just in the unwinding events, like in Eq. 114(1. were performed, for example, in [Durrer 
& Zhou, 1995]. 13 



III. CURRENTS ALONG STRINGS 



In the past few years it has become clear that topological defects, and in particular strings, will be endowed with a 
considerably richer structure than previously envisaged. In generic grand unified models the Higgs field, responsible 
for the existence of cosmic strings, will have interactions with other fundamental fields. This should not surprise us, 
for well understood low energy particle theories include field interactions in order to account for the well measured 
masses of light fermions, like the familiar electron, and for the masses of gauge bosons W and Z discovered at CERN in 
the eighties. Thus, when one of these fundamental (electromagnetically charged) fields present in the model condenses 
in the interior space of the string, there will appear electric currents flowing along the string core. 

Even though these strings are the most attractive ones, the fact of them having electromagnetic properties is not 
actually fundamental for understanding the dynamics of circular string loops. In fact, while in the uncharged and non 
current-carrying case symmetry arguments do not allow us to distinguish the existence of rigid rotations around the 
loop axis, the very existence of a small current breaks this symmetry, marking a definite direction, which allows the 
whole loop configuration to rotate. This can also be viewed as the existence of spinning particle-like solutions trapped 
inside the core. The stationary loop solutions where the string tension gets balanced by the angular momentum of 
the charges is what Davis and Shellard [1988] dubbed vortons. 

Vorton configurations do not radiate classically. Because they have loop shapes, implying periodic boundary 
conditions on the charged fields, it is not surprising that these configurations are quantized. At large distances these 
vortons look like point masses with quantized electric charge (actually they can have more than a hundred times the 
electron charge) and angular momentum. They are very much like particles, hence their name. They are however 
very peculiar, for their characteristic size is of order of their charge number (around a hundred) times their thickness, 
which is essentially some fourteen orders of magnitude smaller than the classical electron radius. Also, their mass is 
often of the order of the energies of grand unification, and hence vortons would be some twenty orders of magnitude 
heavier than the electron. 



In fact, in the action the coupling constant A of the 'Mexican hat' potential is interpreted as the Lagrange multiplier. 
Simulations of the collapse of 'exotic' textures can be found at http: //camelot .mssm. edu/~ats/texture . html 
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But why should strings become conducting in the first place? The physics inside the core of the string differs 
somewhat from outside of it. In particular the existence of interactions among the Higgs field forming the string 
and other fundamental fields, like that of charged fermions, would make the latter loose their masses inside the core. 
Then, only small energies would be required to produce pairs of trapped fermions and, being effectively massless 
inside the string core, they would propagate at the speed of light. These zero energy fermionic states, also called zero 
modes, endow the string with currents and in the case of closed loops they provide the mechanical angular momentum 
support necessary for stabilizing the contracting loop against collapse. 



A. Goto— Nambu Strings 

Our aim now would be to introduce extra fields into the problem and see what new features arise. We would expect 
to find -among other novelties- currents flowing along the defect cores, as advertised before. However, doing this in 
detail would unfortunately take us too much away from the main topic of these notes, and we just refer the reader 
to some recent work [Lemperiere & Shellard 2002] (see also H3,E1) anc ^ *° ^ ne recent review in ^cji as weu as to 
the other references given in the introduction, for a detailed treatment. Here below, we just give some few additional 
features of cosmic string field theory. 

The simple Lagrangian we saw in previous sections was a good approximation for ideal structureless strings, known 
under the name of Goto-Nambu strings [Goto, 1971; Nambu, 1970]. Additional fields coupled with the string-forming 
Higgs field often lead to interesting effects in the form of generalized currents flowing along the string core. But before 
taking into full consideration the internal structure of strings (given in 40] ) it is appropriate to start by setting the 
scene with the simple Abelian Higgs model (which describes scalar electrodynamics). This is a prototype of gauge 
field theory with spontaneous symmetry breaking G = U(l) — * {1}. The Lagrangian reads [Higgs, 1964] 

£ H = -\[D^][D^Y \{F$? - ^(|$| 2 - v 2 ) 2 , (16) 

with gauge covariant derivative = d^ + iqA^, antisymmetric tensor = V^A^ — V„A^ for the gauge vector 
field A„, and complex scalar field $ = |$|e m with gauge coupling q. 

The first solutions for this theory were found by Nielsen & Olesen [1973]. A couple of relevant properties are 
noteworthy: 

• the mass per unit length for the string is \i = U ~ rj 2 . For GUT local strings this gives \i ~ 10 22 g/cm, while 
one finds /i ~ rj 2 ln(r / 'm~ l ) — > oo if strings are global, due to the absence of compensating gauge fields. This 
divergence is in general not an issue, because global strings only in few instances are isolated; in a string network, 
a natural cutoff is the distance to the neighboring string. 

1/2 

• There are essentially two characteristic mass scales (or inverse length scales) in the problem: m s ~ r\ and 

m v ~ qr], corresponding to the inverse of the Compton wavelengths of the scalar (Higgs) and vector (Ajf ) 
particles, respectively. 

• There exists a sort of screening of the energy, called 'Higgs screening', implying a finite energy configuration, 
thanks to the way in which the vector field behaves far from the string core: Ag — ► (l/qr)da/d6 , for r — > oo. 

After a closed path around the vortex one has $(27r) = $(0), which implies that the winding phase a should 
be an integer times the cylindrical angle 0, namely a = nd. This integer n is dubbed the 'winding number'. In 
turn, from this fact it follows that there exists a tube of quantized 'magnetic' flux, given by 

/ r r„ 1 f 27T da ,„ 27m . , 

$ B = J> A.dl = - / —d6 = (17) 

J Q Jo d0 q 

In the string there is a sort of competing effect between the fields: the gauge field acts in a repulsive manner; the 
flux doesn't like to be confined to the core and B lines repel each other. On the other hand, the scalar field behaves 
in an attractive way; it tries to minimize the area where V($>) ^ 0, that is, where the field departs from the true 
vacuum. 

Finally, we can mention a few condensed-matter 'cousins' of Goto-Nambu strings: flux tubes in superconductors 
[Abrikosov, 1957] for the nonrelativistic version of gauge strings ($ corresponds to the Cooper pair wave function). 
Also, vortices in superfluids, for the nonrelativistic version of global strings (<& corresponds to the Bose condensate 
wave function). Moreover, the only two relevant scales of the problem we mentioned above are the Higgs mass m s 
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FIG. 6: Higgs field and energy profiles for Goto-Nambu cosmic strings. The left panel shows the amplitude of the Higgs field 
around the string. The field vanishes at the origin (the false vacuum) and attains its asymptotic value (normalized to unity 
in the figure) far away from the origin. The phase of the scalar field (changing from to 2ir) is shown by the shading of the 
surface. In the right panel we show the energy density of the configuration. The maximum value is reached at the origin, 
exactly where the Higgs is placed in the false vacuum. [Hindmarsh & Kibble, 1995]. 



and the gauge vector mass m v . Their inverse give an idea of the characteristic scales on which the fields acquire their 
asymptotic solutions far away from the string 'location'. In fact, the relevant core widths of the string are given by 
TOjT 1 and m~ 1 . It is the comparison of these scales that draws the dividing line between two qualitatively different 
types of solutions. If we define the parameter [3 — (m s /m v ) 2 , superconductivity theory says that [3 < 1 corresponds to 
Type I behavior while (3 > 1 corresponds to Type II. For us, (3 < 1 implies that the characteristic scale for the vector 
field is smaller than that for the Higgs field and so magnetic field B flux lines arc well confined in the core; eventually, 
an n-vortex string with high winding number n stays stable. On the contrary, (3 > 1 says that the characteristic scale 
for the vector field exceeds that for the scalar field and thus B flux lines are not confined; the n-vortex string will 
eventually split into n vortices of flux 2n/q. In summary: 

a _ i m s x2 / < 1 n— vortex stable (B flux lines confined in core) — Type I , 1s s 
m v \ > 1 Unstable : splitting into n vortices of flux 2ir/q — Type II ^ ' 



IV. STRUCTURE FORMATION FROM DEFECTS 



A. Cosmic strings 



In this section we will provide just a quick description of the remarkable cosmological features of cosmic strings. 
Many of the proposed observational tests for the existence of cosmic strings are based on their gravitational inter- 
actions. In fact, the gravitational field around a straight static string is very unusual [Vilenkin, 1981]. As is well 
known, the Newtonian limit of Einstein field equations with source term given by T{? = diag( i o, — pi, — p 2 , — P3) in 
terms of the Newtonian potential <I> is given by V 2 $ = 4irG(p + pi + P2 + P3), just a statement of the well known 
fact that pressure terms also contribute to the 'gravitational mass'. For an infinite string in the z-direction one has 
P3 = —p, i.e., strings possess a large relativistic tension (negative pressure). Moreover, averaging on the string core 
results in vanishing pressures for the x and y directions yielding V 2< f> = for the Poisson equation. This indicates 
that space is flat outside of an infinite straight cosmic string and therefore test particles in its vicinity should not feel 
any gravitational attraction. 

In fact, a full general relativistic analysis confirms this and test particles in the space around the string feel no 
Newtonian attraction; however there exists something unusual, a sort of wedge missing from the space surrounding 
the string and called the 'deficit angle', usually noted A, that makes the topology of space around the string that of 
a cone. To see this, consider the metric of a source with energy-momentum tensor [Vilenkin 1981, Gott 1985] 

T; = ^)%)diag(/x,0,0,T) . (19) 

In the case with T = p, (a rather simple equation of state) this is the effective energy-momentum tensor of an 
unperturbed string with string tension fi as seen from distances much larger than the thickness of the string (a 
Goto-Nambu string). However, real strings develop small-scale structure and are therefore not well described by 
the Goto-Nambu action. When perturbations are taken into account T and /j, are no longer equal and can only 
be interpreted as effective quantities for an observer who cannot resolve the perturbations along its length. And 
in this case we are left without an effective equation of state. Carter [1990] has proposed that these 'noisy' strings 
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should be such that both its speeds of propagation of perturbations coincide. Namely, the transverse (wiggle) speed 
c T = (T/^i) 1 / 2 for extrinsic perturbations should be equal to the longitudinal (woggle) speed cl = (—dT j dii) 1 ! 2 for 
sound-type perturbations. This requirement yields the new equation of state 

M T = /x 2 (20) 

and, when this is satisfied, it describes the energy-momentum tensor of a wiggly string as seen by an observer who 
cannot resolve the wiggles or other irregularities along the string [Carter 1990, Vilenkin 1990]. 




FIG. 7: Cosmic strings affect surrounding spacetime by removing a small angular wedge, creating a conelike geometry (left). 
Space remains Bat everywhere, but a circular path around the string encompasses slightly less than 360 degrees. The deficit 
angle is tiny, about 10~ radian. To an observer, the presence of a cosmic string would be betrayed by its effect on the trajectory 
of passing light rays, which are deflected by an amount equal to the deficit angle. The resultant gravitational lensing reveals 
itself in the doubling of images of objects behind the string (right panel). 

The gravitational field around the cosmic string [neglecting terms of order (G/i) 2 ] is found by solving the linearized 
Einstein equations with the above T". One gets 

/ioo = ^33 = 4G(^ - T) Hr/ro), (21) 

h 11 = h 22 =4G(fi + T)ln(r/r ), (22) 

where /i M „ = g M „ — r/^ is the metric perturbation, the radial distance from the string is r = (x 2 + y 2 ) 1 ^ 2 , and r$ is a 
constant of integration. 

For an ideal, straight, unperturbed string, the tension and mass per unit length are T = /j, = fiQ and one gets 

h 00 = h 33 = 0, /in = h 2 2 = 8Ghq ln(r/r ). (23) 

By a coordinate transformation one can bring this metric to a locally flat form 

ds 2 = dt 2 - dz 2 - dr 2 - (1 - 8Gu. Q )r 2 d(f> 2 , (24) 

which describes a conical and flat (Euclidean) space with a wedge of angular size A = 8irGn Q (the deficit angle) 
removed from the plane and with the two faces of the wedge identified. 



1. Wakes and gravitational lensing 

We saw above that test particles 14 at rest in the spacetime of the straight string experience no gravitational force, 
but if the string moves the situation radically changes. Two particles initially at rest while the string is far away, will 



If one takes into account the own gravitational field of the particle living in the spacetime around a cosmic string, then the situation 
changes. In fact, the presence of the conical 'singularity' introduced by the string distorts the particle's own gravitational field and 
results in the existence of a weak attractive force proportional to G 2 fim 2 /r 2 , where m is the particle's mass [Linet, 1986]. 
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suddenly begin moving towards each other after the string has passed between them. Their head-on velocities will 
be proportional to A or, more precisely, the particles will get a boost v = 47rG/iow s 7 in the direction of the surface 
swept out by the string. Here, 7 = (1 — v 2 )^ 1 ^ 2 is the Lorentz factor and v s the velocity of the moving string. Hence, 
the moving string will built up a wake of particles behind it that may eventually form the 'seed' for accreting more 
matter into sheet-like structures [Silk & Vilenkin 1984]. 




FIG. 8: By deflecting the trajectory of ordinary matter, strings offer an interesting means of forming large-scale structure. A 
string sweeping through a distribution of interstellar dust will draw particles together in its wake, giving them lateral velocities 
of a few kilometers per second. The trail of the moving string will become a planar region of high-density matter, which, 
after gravitational collapse, could turn into thin, sheetlike distributions of galaxies [fmage courtesy of Pedro Avelino and Paul 
Shellard]. 

Also, the peculiar topology around the string makes it act as a cylindric gravitational lens that may produce 
double images of distant light sources, e.g., quasars. The angle between the two images produced by a typical GUT 
string would be oc Gfi and of order of a few arcseconds, independent of the impact parameter and with no relative 
magnification between the images [see Cowie & Hu, 1987, for a first observational attempt]. Surprisingly enough, a 
recent detection of an extragalactic double source with the appropriate characteristics (few arcseconds angle, lack of 
excess light in between the images, etc), lead Sazhin et al. [2003] to propose the serendipitous discovery of a cosmic 
string lens event. Their data suggest that both images belong to early-type giant elliptical galaxies with redshift 0.46 
(or some 1900 Mpc away, for a reduced Hubble constant h — 0.65) meaning that these galaxies (if indeed there is such 
chance projection) are some 20 Kpc away from each other. On the contrary, if this pair is caused by the splitting due 
to an intervening cosmic string, the energy-scale of the symmetry-breaking transition giving bith to such string can 
be computed, and it turns out to be of a typical GUT scale. Doubtless, more independent observations are needed to 
confirm this interesting case. 

Turning back now to the peculiar matter structures generated by moving strings, the situation described above 
gets even more interesting when we allow the string to have small-scale structure, which we called wiggles, as in 
fact simulations indicate. Wiggles not only modify the string's effective mass per unit length, [a, but also built up 
a Newtonian attractive term in the velocity boost inflicted on nearby test particles. To see this, let us consider the 
formation of a wake behind a moving wiggly string. Assuming the string moves along the a;-axis, we can describe the 
situation in the rest frame of the string. In this frame, it is the particles that move, and these flow past the string 
with a velocity v s in the opposite direction. Using conformally Minkowskian coordinates we can express the relevant 
components of the metric as 

ds 2 = (1 + h o)[dt 2 - {dx 2 + dy 2 )], (25) 
where the missing wedge is reproduced by identifying the half-lines y — ±AnGfix, x > 0. The linearized geodesic 
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equations in this metric can be written as 

2x = -(l-i 2 -y 2 )d x h 00 , (26) 

2y = -{l-x 2 -y 2 )d y hao, (27) 

where over-dots denote derivatives with respect to t. Working to first order in G/z, the second of these equations can 
be integrated over the unperturbed trajectory x — v s t, y = yo. Transforming back to the frame in which the string 
has a velocity v s yields the result for the velocity impulse in the y-direction after the string has passed [Vachaspati 
& Vilenkin, 1991; Vollick, 1992] 

v = _ ^G0m-T) _ 

The second term is the velocity impulse due to the conical deficit angle we saw above. This term will dominate 
for large string velocities, case in which big planar wakes are predicted. In this case, the string wiggles will produce 
inhomogcncities in the wake and may easy the fragmentation of the structure. The 'top-down' scenario of structure 
formation thus follows naturally in a universe with fast-moving strings. On the contrary, for small velocities, it is the 
first term that dominates over the deflection of particles. The origin of this term can be easily understood [Vilenkin & 
Shcllard, 2000]. From Eqn. Q2ip. the gravitational force on a non-relativistic particle of mass m is F ~ mG(p, — T)/r. 
A particle with an impact parameter r is exposed to this force for a time At ~ r/v s and the resulting velocity is 
v ~ (F/m)At ~ G(jjl - T)/v s . 

B. Textures 

During the radiation era, and when the correlation length is already growing with the Hubble radius, the texture 
field has energy density ptexture ~ (V0) 2 ~ rj 2 /H~ 2 , and remains a fixed fraction of the total density p c ~ t~ 2 yielding 
^texture ~ Gij 2 . This is the scaling behavior for textures and thus we do not need to worry about textures dominating 
the universe. 

But as we already mentioned, textures are unstable to collapse, and this collapse generates perturbations in the 
metric of spacetime that eventually lead to large scale structure formation. These perturbations in turn will affect 
the photon geodesies leading to CMB anisotropies, the clearest possible signature to probe the existence of these 
exotic objects being the appearance of hot and cold spots in the microwave maps. Due to their scaling behavior, 
the density fluctuations induced by textures on any scale at horizon crossing are given by [5p/ p)h ~ Grj 2 . CMB 
temperature anisotropies will be of the same amplitude. Numerically- simulated maps, with patterns smoothed over 
10° angular scales, by Bennett & Rhie [1993] yield, upon normalization to the COBE-DMR data, a dimensionlcss 
value Grj 2 ~ 10~ 6 , in good agreement with a GUT phase transition energy scale. It is fair to say, however, that the 
texture scenario is having problems in matching current data on smaller scales [see, e.g., Durrer, 2000]. 

C. Defects as dark energy 

There is recent mounting evidence that our current universe is being dominated by a unexpectedly large amount of 
dark energy [e.g., Riess et al, 1998; Pcrlmutter et al., 1999]. Recent observations with type la supernovae, together 
with other astrophysical tests, suggest that more than 65 percent of the critical energy density is made up by some 
yet unknown energy component. 

Cosmic defects can also be seen as a novel form of dark energy. For example, a tangled web of cosmic strings with 
fixed mass per unit length, which self-intersects without having reconnection. Non intercommuting strings means no 
production of loops, and therefore the main channel for loosing energy is not active. The model proposed in [Vilenkin, 
1984] has the mean mass density in strings scaling as /^strings K (^ a M) _1 instead of ^strings oc rj 2 t~ 2 as we saw above. 
From this, one has Pstrings^ 2 oc t 1 ! 2 in the radiation-dominated era and Pstrings^ 2 oc i 1 / 3 during matter domination, 
which means that the energy in strings firings grows with time and, after a certain t s , strings would dominate the 
universe. With t s falling in the matter-domination era, we have p s trings / p(t s ) ~ (*eq/*«) (^Aeq) 1 '^/* ~ 1, with 
the background p oc 1/Gt 2 . In the case < z < 3, roughly t s ~ 10 17 sec, with t cq ~ 10 11 sec. and ~ mp/r/ 2 , we 
get Gp ~ 10~ 20 for the characteristic energy scale of these non-intercommuting strings. 

After t s the Friedmann's equation can be cast as (a/a) 2 oc G/9 str ings oc G/ta(t), which implies that the scale factor 
goes as a(t) oc t and then /^strings 1/t 2 . Now, recalling the local energy conservation law p = — 3(p + p)a/a, and 
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applying it for a dark u x" component, w x = p x /px, we get p x cx a~ 3 ( 1+W:c \ If this dark component is made up by 
strings, one then deduces that it should be w x — —1/3. Of course, this gives a — for the scale factor, so it cannot 
explain the recent acceleration phase. It nevertheless goes in the right direction. 

Similar arguments have been studied for other defects, like textures [Davis, 1987] and can also be devised for domain 
walls [Zel'dovich et al., 1974; Battye et al, 1999], in this latter case yielding w x = —2/3 which points closer to the 
observational "equation of state" currently selected by the analysis of the different astrophysical surveys. For these 
and other reasons, with the words of the recent authoritative review by Peebles & Ratra [2002], the class of cosmic 
defect models is worth bearing in mind. 

V. CMB SIGNATURES FROM DEFECTS 

If cosmic defects have really formed in the early universe and some of them are still within our present horizon 
today, the anisotropics in the CMB they produce would have a characteristic signature. Strings, for example, would 
imprint the background radiation in a very particular way due to the Doppler shift that the background radiation 
suffers when a string intersects the line of sight. The conical topology of space around the string will produce a 
differential redshift of photons passing on different sides of it, resulting in step-like discontinuities in the effective 
CMB temperature, given by ^f- « %ttG^lv s ^ with, as before, 7 = (1 — v 2 )^ 1 / 2 the Lorentz factor and v s the velocity 
of the moving string. This 'stringy' signature was first studied by Kaiser & Stebbins [1984] and Gott [1985] (see 
Figure EJ). 




FIG. 9: The Kaiser-Stebbins effect for cosmic strings. A string network evolves into a self-similar scaling regime, perturbing 
matter and radiation during its evolution. The effect on the CMB after recombination leads to distinct steplike discontinuities 
on small angular scales that were Erst studied by Kaiser & Stebbins [1984]. The left panel shows a simulated patch of the sky 
that fits in one of the pixels of the COBE experiment. Hence, higher resolution observatories are needed in order to detect 
strings. The right panel shows a patch on the CMB sky of order 20' across. However, recent studies indicate that this clean 
tell-tale signal gets obscured at subdegree angular scales due to the temperature fluctuations generated before recombination. 
[Magueijo & Ferreira 1997]. 

Anisotropics of the CMB are directly related to the origin of structure in the universe. Galaxies and clusters of 
galaxies eventually formed by gravitational instability from primordial density fluctuations, and these same fluctua- 
tions left their imprint on the CMB. Recent balloon [de Bernardis, et ah, 2000; Hanany, et al, 2000], ground-based 
interferometer [Halverson, et al., 2001] and satellite [Bennett, et al. 2003] experiments have produced reliable estimates 
of the power spectrum of the CMB temperature anisotropies. While they helped eliminate certain candidate theories 
for the primary source of cosmic perturbations, the power spectrum data is still compatible with the theoretical es- 
timates of a relatively large variety of models, such as ACDM, quintessence models or some hybrid models including 
cosmic defects. 

There are two main classes of models of structure formation -passive and active models. In passive models, density 
inhomogeneities are set as initial conditions at some early time, and while they subsequently evolve as described 
by Einstein-Boltzmann equations, no additional perturbations are seeded. On the other hand, in active models the 
sources of density perturbations are time-dependent. 

All specific realizations of passive models are based on the idea of inflation. In simplest inflationary models it is 
assumed that there exists a weakly coupled scalar field <fi, called the inflaton, which "drives" the (quasi) exponential 
expansion of the universe. The quantum fluctuations of (f> are stretched by the expansion to scales beyond the horizon, 
thus "freezing" their amplitude. Inflation is followed by a period of thermalization, during which standard forms of 
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matter and energy are formed. Because of the spatial variations of <f> introduced by quantum fluctuations, thermaliza- 
tion occurs at slightly different times in different parts of the universe. Such fluctuations in the thcrmalization time 
give rise to density fluctuations. Because of their quantum nature and because of the fact that initial perturbations 
are assumed to be in the vacuum state and hence well described by a Gaussian distribution, perturbations produced 
during inflation are expected to follow Gaussian statistics to a high degree [Gangui, Lucchin, Matarrese & Mollerach, 
1994], or either be products of Gaussian random variables. This is a fairly general prediction that is being tested 
currently with data from WMAP and will be tested more thoroughly in the future with Planck. 15 

Active models of structure formation are motivated by cosmic topological defects with the most promising candidates 
being cosmic strings. As we saw in previous sections, it is widely believed that the universe underwent a series of 
phase transitions as it cooled down due to the expansion. If our ideas about grand unification are correct, then some 
cosmic defects should have formed during phase transitions in the early universe. Once formed, cosmic strings could 
survive long enough to seed density perturbations. Defect models possess the attractive feature that they have no 
parameter freedom, as all the necessary information is in principle contained in the underlying particle physics model. 
Generically, perturbations produced by active models are not expected to be Gaussian distributed [Gangui, Pogosian 
& Winitzki, 2001a]. 

A. CMB power spectrum from strings 

The narrow main peak and the presence of the second and the third peaks in the CMB angular power spectrum, 
as measured by BOOMERANG, MAXIMA, DASI and WMAP [de Bernardis, et al, 2000; Hanany, et al., 2000; 
Halverson, et al., 2001; Page, et al., 2003], is an evidence for coherent oscillations of the photon-baryon fluid at the 
beginning of the decoupling epoch [see, e.g., Gangui, 2001]. While such coherence is a property of all passive model, 
realistic cosmic string models produce highly incoherent perturbations that result in a much broader main peak. This 
excludes cosmic strings as the primary source of density fluctuations unless new physics is postulated, e.g. models with 
a varying speed of light [Avelino & Martins, 2000]. In addition to purely active or passive models, it has been recently 
suggested that perturbations could be seeded by some combination of the two mechanisms. For example, cosmic 
strings could have formed just before the end of inflation and partially contributed to seeding density fluctuations. It 
has been shown [Contaldi, et al., 1999; Battye & Weller, 2000; Bouchet, et al, 2001] that, although not compelling, 
such hybrid models can be rather successful in fitting the CMB power spectrum data. 

Calculating CMB anisotropies sourced by topological defects is a rather difficult task. In inflationary scenario the 
entire information about the seeds is contained in the initial conditions for the perturbations in the metric. In the case 
of cosmic defects, perturbations are continuously seeded over the period of time from the phase transition that had 
produced them until today. The exact determination of the resulting anisotropy requires, in principle, the knowledge 
of the energy-momentum tensor [or, if only two point functions are being calculated, the unequal time correlators, 
Pen, Seljak, & Turok, 1997] of the defect network and the products of its decay at all times. This information is simply 
not available. Instead, a number of clever simplifications, based on the expected properties of the defect networks {e.g. 
scaling), are used to calculate the source. The latest data from BOOMERANG, MAXIMA and WMAP experiments 
clearly disagree with the predictions of these simple models of defects [Durrer, Gangui & Sakellariadou, 1996]. 

The shape of the CMB angular power spectrum is determined by three main factors: the geometry of the universe, 
coherence and causality. The curvature of the universe directly affects the paths of light rays coming to us from the 
surface of last scattering. In a closed universe, because of the lensing effect induced by the positive curvature, the 
same physical distances between points on the sky would correspond to larger angular scales. As a result, the peak 
structure in the CMB angular power spectrum would shift to larger angular scales or, equivalently, to smaller values 
of the multipoles Ps. 

The prediction of the cosmic string model of [Pogosian & Vachaspati, 1999] for f^totai = 1.3 is shown in Figure ITU1 
As can be seen, the main peak in the angular power spectrum can be matched by choosing a reasonable value for 
^totai- However, even with the main peak in the right place the agreement with the data is far from satisfactory. The 
peak is significantly wider than that in the data and there is no sign of a rise in power at multipole roughly I ss 400 as 
the actual data from WMAP suggests. The sharpness and the height of the main peak in the angular spectrum can 
be enhanced by including the effects of gravitational radiation [Contaldi, Hindmarsh & Magueijo, 1999] and wiggles 
[Pogosian & Vachaspati, 1999]. More precise high-resolution numerical simulations of string networks in realistic 
cosmologies with a large contribution from are needed to determine the exact amount of small-scale structure on 



Useful CMB resources can be found at http://www.mpa-garching.rnpg.de/-banday/CMB.html 
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the strings and the nature of the products of their decay [Landriau & Shellard, 2002]. It is, however, unlikely that 
including these effects alone would result in a sufficiently narrow main peak and some presence of a second peak. 



FIG. 10: The CMB power spectrum produced by the wiggly string model of [Pogosian & Vachaspati, 1999] in a closed universe 
with fitotai = 1.3, fibaryon = 0.05, ficDM = 0.35, Qa = 0.9, and Ho = 65 km s _1 Mpc _1 [Pogosian, 2001]. 

This brings us to the issues of causality and coherence and how the random nature of the string networks comes into 
the calculation of the anisotropy spectrum. Both experimental and theoretical results for the CMB power spectra 
involve calculations of averages. When estimating the correlations of the observed temperature anisotropies, it is 
usual to compute the average over all available patches on the sky. When calculating the predictions of their models, 
theorists find the average over the ensemble of possible outcomes provided by the model. 

In inflationary models, as in all passive models, only the initial conditions for the perturbations are random. The 
subsequent evolution is the same for all members of the ensemble. For wavelengths higher than the Hubble radius, the 
linear evolution equations for the Fourier components of such perturbations have a growing and a decaying solution. 
The modes corresponding to smaller wavelengths have only oscillating solutions. As a consequence, prior to entering 
the horizon, each mode undergoes a period of phase "squeezing" which leaves it in a highly coherent state by the time 
it starts to oscillate. Coherence here means that all members of the ensemble, corresponding to the same Fourier mode, 
have the same temporal phase. So even though there is randomness involved, as one has to draw random amplitudes 
for the oscillations of a given mode, the time behavior of different members of the ensemble is highly correlated. The 
total spectrum is the ensemble-averaged superposition of all Fourier modes, and the predicted coherence results in 
an interference pattern seen in the angular power spectrum as the well-known acoustic peaks. 

In contrast, the evolution of the string network is highly non-linear. Cosmic strings are expected to move at 
relativistic speeds, self-intersect and reconnect in a chaotic fashion. The consequence of this behavior is that the 
unequal time correlators of the string energy-momentum vanish for time differences larger than a certain coherence 
time (t c in Figure ITT|) . Members of the ensemble corresponding to a given mode of perturbations will have random 
temporal phases with the "dice" thrown on average once in each coherence time. The coherence time of a realistic 
string network is rather short. As a result, the interference pattern in the angular power spectrum is completely 
washed out. 

Causality manifests itself, first of all, through the initial conditions for the string sources, the perturbations in 
the metric and the densities of different particle species. If one assumes that the defects are formed by a causal 
mechanism in an otherwise smooth universe then the correct initial condition are obtained by setting the components 
of the stress-energy pseudo-tensor t^ u to zero [Veeraraghavan & Stebbins, 1990; Pen, Spergel & Turok, 1994]. These 
are the same as the isocurvature initial conditions [Hu, Spergel & White, 1997]. A generic prediction of isocurvature 
models (assuming perfect coherence) is that the first acoustic peak is almost completely hidden. The main peak is 
then the second acoustic peak and in flat geometries it appears at I 300 — 400. This is due to the fact that after 
entering the horizon a given Fourier mode of the source perturbation requires time to induce perturbations in the 
photon density. Causality also implies that no superhorizon correlations in the string energy density are allowed. The 
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correlation length of a "realistic" string network is normally between 0.1 and 0.4 of the horizon size. 

An interesting study was performed by Magueijo, Albrecht, Ferreira & Coulson [1996], where they constructed a 
toy model of defects with two parameters: the coherence length and the coherence time. The coherence length was 
taken to be the scale at which the energy density power spectrum of the strings turns from a power law decay for 
large values of k into a white noise at low k. This is essentially the scale corresponding to the correlation length of the 
string network. The coherence time was defined in the sense described in the beginning of this section, in particular, 
as the time difference needed for the unequal time correlators to vanish. Their study showed (see Figure lllfl that by 




FIG. 11: The predictions of the toy model of Magueijo, et al. [1996] for different values of parameters x c , the coherence length, 
and T c , the coherence time. x c oc n/\ c {n), where n is the conformal time and \ c {n) is the correlation length of the network at 
time n. One can obtain oscillations in the CMB power spectrum by fixing either one of the parameters and varying the other. 

accepting any value for one of the parameters and varying the other (within the constraints imposed by causality) 
one could reproduce the oscillations in the CMB power spectrum. Unfortunately for cosmic strings, at least as we 
know them today, they fall into the parameter range corresponding to the upper right corner in Figure ITTI 

In order to get a better fit to present-day observations, cosmic strings must either be more coherent or they have to 
be stretched over larger distances, which is another way of making them more coherent. To understand this imagine 
that there was just one long straight string stretching across the universe and moving with some given velocity. The 
evolution of this string would be linear and the induced perturbations in the photon density would be coherent. By 
increasing the correlation length of the string network we would move closer to this limiting case of just one long 
straight string and so the coherence would be enhanced. 

The question of whether or not defects can produce a pattern of the CMB power spectrum similar to, and including 
the acoustic peaks of, that produced by the adiabatic inflationary models was repeatedly addressed in the literature 
[Contaldi, Hindmarsh & Magueijo 1999; Magueijo, et al. 1996; Liddle, 1995; Turok, 1996; Avelino & Martins, 2000]. In 
particular, it was shown [Magueijo, et al. 1996; Turok, 1996] that one can construct a causal model of active seeds which 
for certain values of parameters can reproduce the oscillations in the CMB spectrum. The main problem today is that 
current realistic models of cosmic strings fall out of the parameter range that is needed to fit the observations. At the 
moment, only the (non-minimal) models with either a varying speed of light or hybrid contribution of strings+inflation 
are the only ones involving topological defects that to some extent can match the observations 16 . One possible way to 
distinguish their predictions from those of inflationary models would be by computing key non-Gaussian statistical 
quantities, such as the CMB bispectrum. 



With the recent first- year WMAP data, specially the cross-correlation between the temperature and polarization maps (TE correlation), 
the situation gets worse for the whole class of causal models. In fact, there exists a clear nontrivial cross-correlation signature for angular 
scales above the size of the acoustic horizon [as shown for example in Kogut, et al. 2003] which is virtually impossible to reproduce with 
any class of non contrived causal models, cosmic defects included. 
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B. CMB bispectrum from active models 

Different cosmological models differ in their predictions for the statistical distribution of the anisotropies beyond 
the power spectrum. Current WMAP data and the future Planck satellite mission will provide high-precision data 
allowing definite estimates of non-Gaussian signals in the CMB. It is therefore important to know precisely which 
are the predictions of all candidate models for the statistical quantities that will be extracted from the new data and 
identify their specific signatures. 

Of the available non-Gaussian statistics, the CMB bispectrum, or the three-point function of Fourier components 
of the temperature anisotropy, has been perhaps the one best studied in the literature [Gangui & Martin, 2000a]. 
There are a few cases where the bispectrum may be deduced analytically from the underlying model. The bispectrum 
can be estimated from simulated CMB sky maps; however, computing a large number of full-sky maps resulting from 
defects is a much more demanding task. Recently a precise numerical code to compute it, not using CMB maps and 
similar to the CMBFAST code 17 for the power spectrum, was developed in [Gangui, Pogosian & Winitzki 2001b]. 
What follows below is an account of this work. 

In a few words, given a suitable model, one can generate a statistical ensemble of realizations of defect matter 
perturbations. We used a modified Boltzmann code based on CMBFAST to compute the effect of these perturbations 
on the CMB and found the bispectrum estimator for a given realization of sources. We then performed statistical 
averaging over the ensemble of realizations to compute the expected CMB bispectrum. (The CMB power spectrum 
was also obtained as a byproduct.) As a first application, we then computed the expected CMB bispectrum from a 
model of simulated string networks first introduced by Albrecht et al. [1997] and further developed in [Pogosian & 
Vachaspati, 1999] and in [Gangui, Pogosian & Winitzki 2001]. 

We assume that, given a model of active perturbations, such as a string simulation, we can calculate the energy- 
momentum tensor X^^x, r) for a particular realization of the sources in a finite spatial volume Vq. Here, x is a 
3-dimensional coordinate and r is the cosmic time. Many simulations are run to obtain an ensemble of random 
realizations of sources with statistical properties appropriate for the given model. The spatial Fourier decomposition 
of T M „ can be written as 

T F (x,r)=^e^(k,r) e lkx , (29) 

k 

where k are discrete. If V is sufficiently large we can approximate the summation by the integral 

^ev(k,rK kx « ^ J d 3 ke^(k,r)e ikx , (30) 

and the corresponding inverse Fourier transform will be 

6^(k,T) = -l/ rf 3 xT F (x,r) e - kx . (31) 

Of course, the final results, such as the CMB power spectrum or bispectrum, do not depend on the choice of Vb- To 
ensure this independence, we shall keep Vb in all expressions where it appears below. 

It is conventional to expand the temperature fluctuations over the basis of spherical harmonics, 

AT/T(n) = ^a ;m y ;m (n), (32) 

lm 

where n is a unit vector. The coefficients a/ m can be decomposed into Fourier modes, 

aim = j§r H) ' 47r / d3k A ' (k) YLi ^' (33) 

Given the sources M „(k, r), the quantities A;(k) are found by solving linearized Einstein-Boltzmann equations and 
integrating along the line of sight, using a code similar to CMBFAST [Seljak & Zaldarriaga, 1996]. This standard 
procedure can be written symbolically as the action of a linear operator Bf v (k) on the source energy-momentum 



http : //physics .nyu. edu/matiasz/CMBFAST/ cmbf ast .html 
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tensor, A/(k) = B^ v (k)Q^{k, t), so the third moment of A;(k) is linearly related to the three-point correlator of 
MJ/ (k, t). Below we consider the quantities A ; (k), corresponding to a set of realizations of active sources, as given. 
The numerical procedure for computing Aj(k) was developed in [Albrecht et al. 1997] and in [Pogosian & Vachaspati, 
1999]. 

The third moment of ai m , namely (ai imi ai 2m2 ai 3m3 ) , can be expressed as 

H y 1+ / 2 +/ 3 (4tt) 3 ^ |d 3 k 1 d 3 k 2 d 3 k 3 ^: mi (k 1 )y i : m2 (k 2 )y^ m3 (k 3 ) (A /l (k 1 ) A i2 (k 2 ) A i3 (k 3 )) . (34) 

A straightforward numerical evaluation of Eq. I|34[) from given sources A; (k) is prohibitively difficult, because 
it involves too many integrations of oscillating functions. However, we shall be able to reduce the computation to 
integrations over scalars [a similar method was employed in Komatsu & Spergel, 2001 and in Wang & Kamionkowski, 
2000]. Due to homogeneity, the 3-point function vanishes unless the triangle constraint is satisfied, 

ki + k 2 + k 3 = 0. (35) 

We may write 

(A h (ki) A h (k 2 ) A, 3 (k 3 )) = 6& (^ + k 2 + k 3 ) P hhh (k 1; k 2 , k 3 ) , (36) 

where the three-point function P/ 1 / 2 ; 3 (ki,k 2 ,k 3 ) is defined only for values of k; that satisfy Eq. (|35|l . Given the 
scalar values k\, fc 2 , k 3 , there is a unique (up to an overall rotation) triplet of directions k^ for which the RHS of 
Eq. H3tjf) does not vanish. The quantity Pi t i 2 i 3 (ki, k 2 , k 3 ) is invariant under an overall rotation of all three vectors k^ 
and therefore may be equivalently represented by a function of scalar values k\ 1 fc 2 , fc 3 , while preserving all angular 
information. Hence, we can rewrite Eq. I|36[) as 

(Aj^ki) A, a (k 2 ) A,,(k3)) =5 (3) (k 1 +k 2 +k 3 )P ili2i3 (fc 1 ,A :2 ,fc 3 ). (37) 
Then, using the simulation volume Vo explicitly, we have 

P hhh (fe, fa, = ^7- (Aix(ki) A /2 (k 2 ) A, 3 (k 3 )) . (38) 

Given an arbitrary direction k x and the magnitudes k± , fc 2 and fc 3 , the directions k 2 and k 3 are specified up to overall 
rotations by the triangle constraint. Therefore, both sides of Eq. (|38|) are functions of scalar ki only. The expression 
on the RHS of Eq. I|38() is evaluated numerically by averaging over different realizations of the sources and over 
permissible directions k^; below we shall give more details of the procedure. 

Substituting Eqs. I|37|l and l|38|l into (|34|l . Fourier transforming the Dirac delta and using the Rayleigh identity, we 
can perform all angular integrations analytically and obtain a compact form for the third moment, 

(ai imi ai 2m2 ai 3m3 ) = 'H\ n j^ m3 J r 2 dr b hhh {r), (39) 

where, denoting the Wigner 3j-symbol by ), we have 



rw _ / (2Zi + l)(2Z 2 + l)(2Z 3 + l) h k h) h h h \ m 
hhh V 4tt ^ J \ m x m 2 m 3 J ' { ' 

and where we have defined the auxiliary quantities 6; 1 / 2 i 3 using spherical Bessel functions ji, 

8 V^ 3 f 

b hhh( r ) = ^3(^)3 J kfdh kjdk 2 kjdk 3 

x 3h(kir)ji 2 {k2r)j h (k 3 r)P h i 2h (k 1 ,k 2 ,k 3 ). (41) 

The volume factor V§ contained in this expression is correct: as shown in the next section, each term A; includes a 

factor Vq~ 2 ^ 3 , while the average quantity Pi 1 i 2 i 3 (k%, k%, kz) oc V^, -3 [cf. Eq. |J3SJ], so that the arbitrary volume Vo of 
the simulation cancels. 

Our proposed numerical procedure therefore consists of computing the RHS of Eq. (|39|l by evaluating the necessary 
integrals. For fixed {hhh}, computation of the quantities bi t i 2 i 3 {r) is a triple integral over scalar ki defined by 
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Eq. I|41[) : it is followed by a fourth scalar integral over r [Eq. (|39() ]. We also need to average over many realizations 
of sources to obtain Pi 1 i 2 i 3 (ki, /C2, ^3). It was not feasible for us to precompute the values Pi 1 i 2 i 3 (ki, ki, k$) on a grid 
before integration because of the large volume of data: for each set {^1^3} the grid must contain ~ 10 3 points for 
each ki. Instead, we precompute A;(k) from one realization of sources and evaluate the RHS of Eq. I|38|) on that 
data as an estimator of Pi 1 i 2 i 3 (ki, &2, fc 3 ), averaging over allowed directions of lq. The result is used for integration 
in Eq. (gTJ. 

Because of isotropy and since the allowed sets of directions k 4 are planar, it is enough to restrict the numerical cal- 
culation to directions ki within a fixed two-dimensional plane. This significantly reduces the amount of computations 
and data storage, since A;(k) only needs to be stored on a two-dimensional grid of k. 

In estimating Pi 1 i 2 i 3 (ki,k2,ks) from Eq. I|38|) . averaging over directions of kj plays a similar role to ensemble 
averaging over source realizations. Therefore if the number of directions is large enough (we used 720 for cosmic 
strings), only a moderate number of different source realizations is needed. The main numerical difficulty is the highly 
oscillating nature of the function 6/ 1 ; 2 ; 3 (r). The calculation of the bispectrum for cosmic strings presented in the next 
Section requires about 20 days of a single-CPU workstation time per realization. 

We note that this method is specific for the bispectrum and cannot be applied to compute higher-order correlations. 
The reason is that higher-order correlations involve configurations of vectors k^ that are not described by scalar values 
ki and not restricted to a plane. For instance, a computation of a 4-point function would involve integration of highly 
oscillating functions over four vectors k^ which is computationally infeasible. 

From Eq. (j39|l we derive the CMB angular bispectrum Ci 1 i 2 i 3 , defined as [Gangui & Martin, 2000b] 

{ai imi ai 2 rn 2 ai 3 m 3 ) = 1 1 ' 2 ^ ) Chi2h ■ (42) 
V mi vn2 rn,3 I 

The presence of the 3j-symbol guarantees that the third moment vanishes unless mi + m 2 + m 3 = and the U indices 
satisfy the triangle rule \U — lj\ < Ik < h + lj- Invariance under spatial inversions of the three-point correlation 
function implies the additional 'selection rule' h+h + h — even, in order for the third moment not to vanish. Finally, 
from this last relation and using standard properties of the 3j-symbols, it follows that the angular bispectrum C; 1 z 2 ; 3 
is left unchanged under any arbitrary permutation of the indices It. 

In what follows we will restrict our calculations to the angular bispectrum C; 1 / 2 / 3 in the 'diagonal' case, i.e. h = 
I2 = I3 = I. This is a representative case and, in fact, the one most frequently considered in the literature. Plots of 
the power spectrum are usually done in terms of 1(1 + 1)C; which, apart from constant factors, is the contribution to 
the mean squared anisotropy of temperature fluctuations per unit logarithmic interval of I. In full analogy with this, 
the relevant quantity to work with in the case of the bispectrum is 



G m = 1(21 + If' 2 ' ' ' \Cm ■ (43) 




For large values of the multipole index I, Gm oc l 3 ^ 2 Cm- Note also what happens with the 3j-symbols appearing 
in the definition of the coefficients W™^™ 2 " 13 : the symbol \m1m2m3) ^ s absent from the definition of C^i 2 i 3 , while in 
Eq. i|4*5|) the symbol (000) i s squared. Hence, there are no remnant oscillations due to the alternating sign of (qoo)- 
However, even more important than the value of Cm itself is the relation between the bispectrum and the cosmic 
variance associated with it. In fact, it is their comparison that tells us about the observability 'in principle' of the 
non-Gaussian signal. The cosmic variance constitutes a theoretical uncertainty for all observable quantities and comes 
about due to the fact of having just one realization of the stochastic process, in our case, the CMB sky [Scaramella 
& Vittorio, 1991]. 

The way to proceed is to employ an estimator Gi t i 2 i s for the bispectrum and compute the variance from it. By 
choosing an unbiased estimator we ensure it satisfies Cj 1 i 2 i 3 = (G^i^). However, this condition does not isolate a 
unique estimator. The proper way to select the best unbiased estimator is to compute the variances of all candidates 
and choose the one with the smallest value. The estimator with this property was computed in [Gangui & Martin, 
2000b] and is 

Chhh = f 1 ^ ^ ) a lim i ai 2 m 2 a hm 3 - (44) 

e — ' \ mi mi mi / 

mi,m 2 ,m 3 \ u / 

The variance of this estimator, assuming a mildly non-Gaussian distribution, can be expressed in terms of the angular 
power spectrum C; as follows 

° 2 c, , ,^ c hGi 2 Ci 3 (l + S h i 2 +5i 2 i 3 +Si 3 i 1 +2Si 1 i 2 Si 2 i 3 ) . (45) 
^i\i 2 i 3 
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The theoretical signal-to-noise ratio for the bispectrum is then given by 

(S/N) hhl3 = \C hl2l Ja G |. (46) 

In turn, for the diagonal case l± = I2 = h = I we have 

(S/N) l = \Cui/a dm \. (47) 

Incorporating all the specifics of the particular experiment, such as sky coverage, angular resolution, etc., will allow 
us to give an estimate of the particular non-Gaussian signature associated with a given active source and, if observable, 
indicate the appropriate range of multipole Z's where it is best to look for it. 

C. CMB bispectrum from strings 

To calculate the sources of perturbations we have used an updated version of the cosmic string model first introduced 
by Albrecht et al. [1997] and further developed in [Pogosian & Vachaspati, 1999], where the wiggly nature of strings 
was taken into account. In these previous works the model was tailored to the computation of the two-point statistics 
(matter and CMB power spectra). When dealing with higher-order statistics, such as the bispectrum, a different 
strategy needs to be employed. 

In the model, the string network is represented by a collection of uncorrelated straight string segments produced 
at some early epoch and moving with random uncorrelated velocities. At every subsequent epoch, a certain fraction 
of the number of segments decays in a way that maintains network scaling. The length of each segment at any time 
is taken to be equal to the correlation length of the network. This and the root mean square velocity of segments 
are computed from the velocity-dependent one-scale model of Martins & Shellard [1996]. The positions of segments 
are drawn from a uniform distribution in space, and their orientations are chosen from a uniform distribution on a 
two-sphere. 

The total energy of the string network in a volume V at any time is E = NfiL, where N is the total number of 
string segments at that time, \l is the mass per unit length, and L is the length of one segment. If L is the correlation 
length of the string network then, according to the one-scale model, the energy density is p = E/V = p/L 2 , where 
V = Vba 3 , the expansion factor a is normalized so that a = 1 today, and Vb is a constant simulation volume. It follows 
that N = V/L 3 = Vo/£ 3 , where I = L/a is the comoving correlation length. In the scaling regime I is approximately 
proportional to the conformal time r and so the number of strings N(t) within the simulation volume Vq falls as t~ 3 . 

To calculate the CMB anisotropy one needs to evolve the string network over at least four orders of magnitude in 
cosmic expansion. Hence, one would have to start with N J> 10 12 string segments in order to have one segment left 
at the present time. Keeping track of such a huge number of segments is numerically infeasible. A way around this 
difficulty was suggested in Ref.Q, where the idea was to consolidate all string segments that decay at the same epoch. 
The number of segments that decay by the (discretized) conformal time Tj is 

N d (Ti)=Vo(n(Ti-i)-n(n)), (48) 

where n(r) = [£(r)]~ 3 is the number density of strings at time r. The energy- momentum tensor in Fourier space, 
9* , of these Nafa) segments is a sum 

N d (Ti) 

m— 1 

where 0^™ is the Fourier transform of the energy-momentum of the m-th segment. If segments are uncorrelated, then 

(e™e™ ) = <w (e™e™) (50) 

and 

<e™e™'e^") = £ m „^ mm ,<e^e™e^>. (51) 

Here the angular brackets (...) denote the ensemble average, which in our case means averaging over many realizations 
of the string network. If we are calculating power spectra, then the relevant quantities are the two-point functions of 
9^„, namely 

Hu Q U) = ( E £ ©PC')' (52) 
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Eq. 1|50|) allows us to write 

(%u Q i P ) = E < m"O = N d { n ){Q'%Q%), (53) 

m— 1 

where 9^ is of the energy-momentum of one of the segments that decay by the time . The last step in Eq. I|53|) 
is possible because the segments are statistically equivalent. Thus, if we only want to reproduce the correct power 
spectra in the limit of a large number of realizations, we can replace the sum in Eq. (|49|) by 

%„ = y/N^e%. (54) 

The total energy-momentum tensor of the network in Fourier space is a sum over the consolidated segments: 

K K 

e^E^^Ev^R©; 1 ,- (55) 

i=l i=l 

So, instead of summing over 2j=i ^d{Ti) ^ 10 12 segments we now sum over only K segments, making K a parameter. 
For the three-point functions we extend the above procedure. Instead of Eqs. (|52|l and l|53|) we now write 

N d (n)N d (Ti) N d (n) N d (n) 

(% v ei P @ U)={ E E E 0™ ™;' ™'')=E (9%e%*%) = N dW{9%9%e%) (56) 

m=l m' = l m" = l m=l 

Therefore, for the purpose of calculation of three-point functions, the sum in Eq. I|49(l should now be replaced by 

% v = [N d (n)] 1/3 % 1 v - (57) 

Both expressions in Eqs. (|54|l and l|57(l . depend on the simulation volume, Vq, contained in the definition of Nd{ri) 
given in Eq. I|48(l . This is to be expected and is consistent with our calculations, since this volume cancels in expressions 
for observable quantities. 

Note also that the simulation model in its present form does not allow computation of CMB sky maps. This is 
because the method of finding the two- and three-point functions as we described involves "consolidated" quantities 
G^j, which do not correspond to the energy-momentum tensor of a real string network. These quantities are auxiliary 
and specially prepared to give the correct two- or three-point functions after ensemble averaging. 

In Fig. HI we show the results for G^ 3 [cf. Eq. g3J)]. It was calculated using the string model with 800 consolidated 
segments in a flat universe with cold dark matter and a cosmological constant. Only the scalar contribution to the 
anisotropy has been included. Vector and tensor contributions are known to be relatively insignificant for local cosmic 
strings and can safely be ignored in this model 0, 19]} 18 ■ The plots are produced using a single realization of the 

string network by averaging over 720 directions of kj. The comparison of G\(^ (or equivalently G^ 3 ) with its cosmic 
variance [cf. Eq. (|45ll ] clearly shows that the bispectrum (as computed from the present cosmic string model) lies 
hidden in the theoretical noise and is therefore undetectable for any given value of I. 

Let us note, however, that in its present stage the string code employed in these computations describes Brownian, 
wiggly long strings in spite of the fact that long strings are very likely not Brownian on the smallest scales, as recent 
field-theory simulations indicate. In addition, the presence of small string loops [Wu, et al., 1998] and gravitational 
radiation into which they decay were not yet included in this model. These are important effects that could, in 
principle, change the above predictions for the string-generated CMB bispectrum on very small angular scales. 

The imprint of cosmic strings on the CMB is a combination of different effects. Prior to the time of recombination 
strings induce density and velocity fluctuations on the surrounding matter. During the period of last scattering these 
fluctuations are imprinted on the CMB through the Sachs- Wolfe effect, namely, temperature fluctuations arise because 
relic photons encounter a gravitational potential with spatially dependent depth. In addition to the Sachs- Wolfe effect, 
moving long strings drag the surrounding plasma and produce velocity fields that cause temperature anisotropies due 
to Doppler shifts. While a string segment by itself is a highly non-Gaussian object, fluctuations induced by string 



The contribution of vector and tensor modes is large in the case of global strings [Turok, Pen & Seljak, 1998; Durrer, Gangui & 
Sakellariadou, 1996]. 
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FIG. 12: The CMB angular bispectrum in the 'diagonal' case (GI/, 3 ) from wiggly cosmic strings in a spatially flat model with 
cosmological parameters S^cdm = 0.3, fibaryon = 0.05, Qa = 0.65, and Hubble constant H — 0.65 kms~ 1 Mpc~ 1 [upper panel] . In 
the lower panel we show the ratio of the signal to theoretical noise \Ciu/o£ | 1//3 for different multipole indices. Normalization 
follows from fitting the power spectrum to the BOOMERANG and MAXIMA data. 



segments before recombination are a superposition of effects of many random strings stirring the primordial plasma. 
These fluctuations are thus expected to be Gaussian as a result of the central limit theorem. 

As the universe becomes transparent, strings continue to leave their imprint on the CMB mainly due to the Kaiser & 
Stebbins [1984] effect. As we mentioned in previous sections, this effect results in line discontinuities in the temperature 
field of photons passing on opposite sides of a moving long string. 19 However, this effect can result in non-Gaussian 
perturbations only on sufficiently small scales. This is because on scales larger than the characteristic inter-string 
separation at the time of the radiation-matter equality, the CMB temperature perturbations result from superposition 
of effects of many strings and are likely to be Gaussian. Avelino et al. [1998] applied several non-Gaussian tests to 
the perturbations seeded by cosmic strings. They found the density field distribution to be close to Gaussian on 
scales larger than 1.5(£1m^ 2 ) -1 Mpc, where £Im is the fraction of cosmological matter density in baryons and CDM 
combined. Scales this small correspond to the multipole index of order I ~ 10 4 . 



D. CMB polarization 



The possibility that the CMB be polarized was first discussed by Martin Rees in 1968, in the context of anisotropic 
universe models. In spite of his optimism, and after many attempts during more than thirty years, including some 
important upper limits [e.g., Keating, et al. 2001; Hedman, et al. 2001, 2002], there has been no positive detection of 
the polarization field until the DASI detection in September 2002 [Leitch et al. 2002; Kovac et al. 2002]. 

Unlike previous experiments, DASI reached the required sensitivity to make a sounding discovery on angular scales 
~ 0?5. Along the same line, WMAP confirmed this detection with a full- sky coverage and polarization data on five 
different frequencies on angular scales bigger than 0?2. Polarization is an important probe both for cosmological 
models and for the more recent history of our nearby Universe. It arises from the interactions of CMB photons with 
free electrons; hence, polarization can only be produced at the last scattering surface (its amplitude depends on the 
duration of the decoupling process) and, unlike temperature fluctuations, it is largely unaffected by variations of the 
gravitational potential after last scattering 20 . Future measurements of polarization will thus provide a clean view of 
the inhomogeneities of the Universe at about 400,000 years after the Bang. 



The extension of the Kaiser-Stebbins effect to polarization will be treated below. In fact, Benabed and Bernardeau [2000] have recently 
considered the generation of a B-type polarization field out of E-type polarization, through gravitational lensing on a cosmic string. 
With the formation of the first stars and quasars, and the subsequent UV radiation emited by these primitive sources, the hydrogen can 
re-ionize. As a consequence, the CMB will scatter again upon ionized matter and will also modify its polarization, albeit on a different 
angular scale. Data from first-year WMAP indicates that reionization did indeed take place somewhere around redshifts z ~ 20 (with 
big "error" bars), which, translated to the elapsed time since the big bang, represents roughly a few hundred millon years. 
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FIG. 13: An electromagnetic linearly polarized wave (in red) oscillates in a given plane (in pink). Reaching an electron (orange 
ball) the wave induces the electron to also oscillate, making it emit radiation (in green). This resulting electromagnetic wave 
is concentrated essentially in the (green ) plane orthogonal to the movement of the electron and it is polarized like the incident 
wave. 

For understanding polarization, a couple of things should be clear. First, the energy of the photons is small compared 
to the mass of the electrons. Then, the CMB frequency does not change, since the electron recoil is negligible. Second, 
the change in the CMB polarization (i.e., the orientation of the oscillating electric field E of the radiation) occurs 
due to a certain transition, called Thomson scattering. The transition probability per unit time is proportional to a 
combination of the old (e a m ) and new (e a out ) directions of polarization in the form |e Q m • e Q out | 2 . In other words, the 
initial direction of polarization will be favored. Third, an oscillating E will push the electron to also oscillate; the 
latter can then be seen as a dipolc (not to be confused with the CMB dipole) , and dipole radiation emits preferentially 
perpendicularly to the direction of oscillation. These 'rules' will help us understand why the CMB should be linearly 
polarized [39| . 

Previous to the recombination epoch, the radiation field is unpolarized. In unpolarizcd light the electric field can 
be decomposed into the two orthogonal directions (along, say, x and z) perpendicular to the line of propagation (y). 
The electric field along e- ln (suppose z is vertical) will make the electron oscillate also vertically. Hence, the dipolar 
radiation will be maximal over the horizontal xy-plane. Analogously, dipole radiation due to the electric field along 
x will be on the yz-plane. If we now look from the side (e.g., from x, on the horizontal plane and perpendicularly 
to the incident direction y) we will see a special kind of scattered radiation. From our position we cannot perceive 
the radiation that the electron oscillating along the x direction would emit, just because this radiation goes to the 
yz-plane, orthogonal to us. Then, it is as if only the vertical component (e- ln ) of the incoming electric field would 
cause the radiation we perceive. From the above rules we know that the highest probability for the polarization of the 
outgoing radiation e Q out will be to be aligned with the incoming one e- m , and therefore it follows that the outgoing 
radiation will be linearly polarized. Now, as both the chosen incoming direction and our position as observers were 
arbitrary, the result will not be modified if we change them. Thomson scattering will convert unpolarized radiation 
into linearly polarized one. 

This however is not the end of the story. To get the total effect we need to consider all possible directions from which 
photons will come to interact with the target electron, and sum them up. We see easily that for an initial isotropic 
radiation distribution the individual contributions will cancel out: just from symmetry arguments, in a spherically 
symmetric configuration no direction is privileged, unlike the case of a net linear polarization which would select one 
particular direction. 

Fortunately, we know the CMB is not exactly isotropic; to the millikelvin precision the dominant mode is dipolar. 
So, what about a CMB dipolar distribution ? Although spatial symmetry does not help us now, a dipole will not 
generate polarization either. Take, for example, the radiation incident onto the electron from the left to be more 
intense than the radiation incident from the right, with average intensities above and below (that's a dipole); it then 
suffices to sum up all contributions to see that no net polarization survives. However, if the CMB has a quadrupolar 
variation in temperature (that it has, first discovered by COBE, to tens of /iK precision) then there will be an excess 
of vertical polarization from left- and right-incident photons (assumed hotter than the mean) with respect to the 
horizontal one from top and bottom light (cooler). From any point of view, orthogonal contributions to the final 
polarization will be different, leaving a net linear polarization in the scattered radiation. 
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FIG. 14: Left panel: non-polarized electromagnetic wave can be decomposed into the sum of two linearly-polarized waves, one 
along the line of sight (in pink), the other along a perpendicular direction (in green). Scattered radiation due to the first wave 
is contained in the plane orthogonal to the line of sight and cannot be detected. Only the second component (in green ) reaches 
the observer and it is polarized as the incident wave. Right panel: when the charged particle receives non polarized waves from 
different directions (in green and pink), it will re-emit the radiation, polarized also along different directions, to the observer. 
If the original radiation is not isotropic (say, the pink arrow is bigger than the green one), then one of the resulting waves (in 
pink) will be slightly more intense than the other (in green), and the observer will perceive a net excess of linear polarization. 



There is one more point to emphasize. Before recombination, ionized matter, electrons and radiation formed a 
single fluid. In it, the inertia was provided by massive nucleons whilst the pressure was that of radiation. And this 
fluid supported sound waves. In fact, the gravitational clumping tendency of the effective mass in the perturbations 
was resisted by the restoring radiation pressure, and therefore gravity-driven acoustic oscillations in both the fluid 
density and local velocity appeared. 

Whereas the acoustic peaks in the temperature anisotropies correspond to the compression and rarefaction maxima 
of the oscillating plasma, the polarization field responds to the local quadrupole moment during the decoupling process. 
But this local quadrupole is mainly due to the Doppler shifts induced by the velocity field of the plasma [Zaldarriaga 
& Harari, 1995]. That is why we know with certainty that polarization shows the uncontaminatcd dynamics of the 
primordial seeds at recombination. 

Within standard recombination models the predicted level of linear polarization on large scales is tiny (see Figure 
1150 : the quadrupole generated in the radiation distribution as the photons travel between successive scatterings is too 
small. Multiple scatterings make the plasma very homogeneous and only wavelengths that are small enough (big fs) 
to produce anisotropies over the (rather short) mean free path of the photons will lead to a significant quadrupole, 
and thus also to polarization. Indeed, if the CMB photons last scattered at z ~ 1100, the SCDM model with h = 1 
predicts no more than 0.05 /xK on scales greater than a few degrees. Hence, measuring polarization at these scales 
represents an experimental challenge. 

However, CMB polarization increases remarkably around the degree-scale in standard models. In fact, for 9 < 1° 
a bump with superimposed acoustic oscillations reaching ~ 5/iK is generically forecasted. On these scales, like for 
the temperature anisotropies, the polarization field shows acoustic oscillations. However, polarization spectra are 
sharper: temperature fluctuations receive contributions from both density (dominant) and velocity perturbations 
and these, being out of phase in their oscillation, partially cancel each other. On the other hand, polarization is 
mainly produced by velocity gradients in the baryon-photon fluid before last scattering, which also explains why 
temperature and polarization peaks are located differently. Moreover, acoustic oscillations depend on the nature 
of the underlying perturbation; hence, we do not expect scalar acoustic sound-waves in the baryon-photon plasma, 
propagating with characteristic adiabatic sound speed cs ~ c/\Z3, close to that of an ideal radiative fluid, to produce 
the same peak-frequency as that produced by gravitational waves, which propagate with the speed of light c (see 

Fignsji. 

The main technical complication with polarization (characterized by a tensor field) is that it is not invariant under 
rotations around a given direction on the sky, unlike the temperature fluctuation that is described by a scalar quantity 
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FIG. 15: Left panel: CMB Polarization for two different models. Red and orange (unlabeled) curves are the angular spectra 
derived for a ACHDM model, both with (red dashed line) and without (red full line) reionization. The temperature anisotropy 
spectrum from scalar perturbations (proportional to [Ce] 1 / 2 , orange curve) is virtually unchanged for both ionization histories. 
The polarization spectrum (cc [C^ ] 1//2 , red curves), although indistinguishable for £ <; 20, dramatically changes for small £'s; 
in this model the Universe is reionized suddenly at low redshift with optical depth r = 0.05 [recall, however, that recent first 
year data from the WMAP satellite indicates that r = 0.17 ± 0.04]. Blue and violet curves represent a SCDM model but with a 
high tensor-mode amplitude, T/S=l at the quadrupole (1 = 2) level, with scale-invariant spectral indices ns = 1 and m = 0. 
Separate scalar (noted C^ S ' ) and tensor (Cy T ' ) contributions to temperature anisotropies are shown (top curves). Scalar modes 
only generate E-type polarization (Cf^), which is smaller than the corresponding red curve of the ACHDM model both due 
to differences in the models (notably A ^ for the red curves) and due to the influence of tensors on the normalization at small 
I. E- and E-type polarization from tensor modes are also shown, respectively C t and C e . Model spectra were computed 
with CMBFAST and are normalized to <5T^ = io = 27.9/iK. Right panel: image of the intensity and polarization of the CMB 
made with the DASI telescope. The small temperature variations of the CMB are shown in false color, with yellow hot and 
red cold. The polarization at each spot in the image is shown by a black line. The length of the line shows the strength of the 
polarization and the orientation of the line indicates the direction in which the radiation is polarized. The size of the white 
spot in the lower left corner approximates the angular resolution of the observations. 

and invariant under such rotations. The level of linear polarization is conveniently expressed in terms of the Stokes 
parameters Q and U. It turns out that there is a clever combination of these parameters that results in scalar quantities 
(in contrast to the above noninvariant tensor description) but with different transformation properties under spatial 
inversions (parity transformations). Then, inspired by classical electromagnetism, any polarization pattern on the 
sky can be separated into 'electric' (scalar, unchanged under parity transformation) and 'magnetic' (pseudo-scalar, 
changes sign under parity) components (E- and B-type polarization, respectively). 

1. CMB polarization from global defects 

One then expands these different components in terms of spherical harmonics, very much like we did for temperature 
anisotropies, getting coefficients a™ for E and B polarizations and, from these, the multipoles Cf' B . The interesting 
thing is that (for symmetry reasons) scalar-density perturbations will not produce any B polarization (a pseudo- 

scalar), that is C e = 0. We see then that an unambiguous detection of some level of B-typc fluctuations will be 
a signature of the existence (and of the amplitude) of a background of gravitational waves ! [Seljak & Zaldarriaga, 
1997] (and, if present, also of rotational modes, like in models with topological defects). 

Linear polarization is a symmetric and traceless 2x2 tensor that requires 2 parameters to fully describe it: Q, U 
Stokes parameters. These depend on the orientation of the coordinate system on the sky. It is convenient to use 
Q + ill and Q — ill as the two independent combinations, which transform under right-handed rotation by an angle 
4> as (Q + ill)' = er 2l<i '(Q + ill) and (Q - ill)' = e 2 '*(Q - ill). These two quantities have spin-weights 2 and -2 
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respectively and can be decomposed into spin ±2 spherical harmonics ±2^r 



(Q + iU)(n) = y j a 2 j m2 Y lm (n) (58) 



I ill 

(Q-iU)(n) = 5>-2,te -2*W«). (59) 
i m 

Spin s spherical harmonics form a complete orthonormal system for each value of s. Important property of spin- 
weighted basis: there exists spin raising and lowering operators 6 and d . By acting twice with a spin lowering 
and raising operator on (Q + iU) and (Q — ill) respectively one obtains quantities of spin 0, which are rotationally 
invariant. These quantities can be treated like the temperature and no ambiguities connected with the orientation 
of coordinate system on the sky will arise. Conversely, by acting with spin lowering and raising operators on usual 
harmonics spin s harmonics can be written explicitly in terms of derivatives of the usual spherical harmonics. Their 
action on ±2^m leads to 

3 2 (Q + t U)(n) = ^(|±|T) ' a 2 , lm Y lm (h) (60) 



lr, 



6 2 (Q~iU){n) = V(£±|1T) ' a. 2 , lm Y lm (n). (61) 



h 



-2]! 



With these definitions the expressions for the expansion coefficients of the two polarization variables become [Seljak 
& Zaldarriaga, 1997] 



0>2,lm — 



0,-llr. 



[l + 2]\ 
[1-2]! 



1/2 



1/2 



dnYC m (h)6 2 {Q + iU){n) (62) 
dnYC m (h)6 2 {Q-iU){h). (63) 



.[Z + 2]! 

Instead of a2,i m , a--i,im it is convenient to introduce their linear electric and magnetic combinations 

1 i 

a-E,lm = —-^{a2,lm + 0,-2,lm) aB,lrn = ^(^.im ~ a -2,lm)- (64) 

These two behave differently under parity transformation: while E remains unchanged B changes the sign, in analogy 
with electric and magnetic fields. 

To characterize the statistics of the CMB perturbations only four power spectra are needed, those for X — T, E, B 
and the cross correlation between T and E. The cross correlation between B and E or B and T vanishes because B 
has the opposite parity of T and E. As usual, the spectra are defined as the rotationally invariant quantities 



Cxi = l^-J ^2{ a X,hn a XJm) C C l = ^2{0*T,lrn a Edm) (65) 



in terms of which on has 



( a X,l'm' a X,lm) — CxiSl'lSm'm (66) 
( a T,l'm' a E,lm) = Cci Sl'lSm'm (67) 
( a *B,l'm' a Edm) = (a* B ,l> m> a T,lm) = 0. (68) 

According to what was said above, one expects some amount of polarization to be present in all possible cosmological 
models. However, symmetry breaking models giving rise to topological defects differ from inflationary models in several 
important aspects, two of which are the relative contributions from scalar, vector and tensor modes and the coherence 
of the seeds sourcing the perturbation equations. In the local cosmic string case one finds that in general scalar modes 
are dominant, if one compares to vector and tensor modes in the usual decomposition of perturbations. The situation 
with global topological defects is radically different and this leads to a very distinctive signature in the polarization 
field. 

Temperature and polarization spectra for various symmetry breaking models were calculated by Seljak, Pen & 
Turok [1997] and are shown in figure 1161 Both electric and magnetic components of polarization are shown for a 
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Waves per radian 



FIG. 16: Left panel: Power spectra of temperature (T), electric type polarization (E) and magnetic type polarization (B) for 
global strings, monopoles, textures and nontopological textures [taken from Seljak. et al., 1997]. The corresponding spectra for 
a standard CDM model with T/S = 1 is also shown for comparison. B polarization turns out to be notably larger for all global 
defects considered if compared to the corresponding predictions of inflationary models on small angular scales. Right panel: 
current and future polarization data by Hivon & Kamionkowski [2002]. Top curve shows the prediction for the polarization 
from primordial inhornogeneities produced by inflation. The large-angle bump in this curve is the enhancement from early star 
formation (reionization) . The lower curves are for inflationary gravitational-wave and gravitational-lensing signals. Recently 
detected DASI data points are shown in red while the rest are expected data points for future experiments with more sensitivity. 

variety of global defects. They also plot for comparison the corresponding spectra in a typical inflationary model, 
namely, the standard CDM model (h — 0.5, = 1, Obaryon = 0.05) but with equal amount of scalars and tensors 
perturbations (noted T/S = 1) which maximizes the amount of B component from inflationary models. In all the 
models they assumed a standard reionization history. The most interesting feature they found is the large magnetic 
mode polarization, with a typical amplitude of ~ lfiK on degree scales [exactly those scales probed by Hedman, et al., 
2001]. For multipoles below £ ~ 100 the contributions from E and B are roughly equal. This differs strongly from the 
inflationary model predictions, where B is much smaller than E on these scales even for the extreme case of T/S ~ 1. 
Inflationary models only generate scalar and tensor modes, while global defects also have a significant contribution 
from vector modes. As we mentioned above, scalar modes only generate E, vector modes predominantly generate B, 
while for tensor modes E and B are comparable with B being somewhat smaller. Together this implies that B can 
be significantly larger in symmetry breaking models than in inflationary models. In figure [TBI we also show the recent 
discovery of a tiny level of polarization by the DASI collaboration together with predictions for future experiments, 
assuming an inflationary origin for the temperature perturbation and polarization signals 21 . 

2. String lensing and CMB polarization 

Recent studies have shown that in realistic models of inflation cosmic string formation seems quite natural in a 
post-inflationary preheating phase [Tkachev et al, 1998, Kasuya & Kawasaki, 1998]. So, even if the gross features 
on CMB maps are produced by a standard (e.g., inflationary) mechanism, the presence of defects, most particularly 
cosmic strings, could eventually leave a distinctive signature. One such feature could be found resorting to CMB 
polarization: the lens effect of a string on the small scale E'-type polarization of the CMB induces a significant 
amount of B-type polarization along the line-of-sight [Zaldarriaga & Seljak, 1998; Benabed & Bernardeau 2000]. This 
is an effect analogous to the Kaiser-Stcbbins effect for temperature maps. 

In the inflationary scenario, scalar density perturbations generate a scalar polarization pattern, given by B-type 
polarization, while tensor modes have the ability to induce both E and B types of polarization. However, tensor 
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modes contribute little on very small angular scales in these models. So, if one considers, say, a standard ACDM 
model, only scalar primary perturbations will be present without defects. But if a few strings are left from a very 
early epoch, by studying the patch of the sky where they are localized, a distinctive signature could come to light. 

In the small angular scale limit, in real space and in terms of the Stokes parameters Q and U one can express the 
E and B fields as follows 

E = A- 1 [(dx 2 -dy 2 )Q + 2dxdyU], (69) 
B = A^^dx 2 -dy 2 )U -2dxdyQ}. (70) 

The polarization vector is parallel transported along the geodesies. The lens affects the polarization by displacing 
the apparent position of the polarized light source. Hence, the observed Stokes parameters Q and U are given in 
terms of the primary (unlcnscd) ones by: Q(a) = Q(a + £) and U(a) = U(a + £). The displacement £ is given by 
the integration of the gravitational potential along the line-of-sights. Of course, here the 'potential' acting as lens is 
the cosmic string whose effect on the polarization field we want to study. 




FIG. 17: Simulations for the B field in the case of a circular loop. The angular size of the figure is 50' x 50'. The resolution is 
5' (left) and 1.2' (right). The discontinuity in the B field is sharper the better the resolution. Weak lensing of CMB photons 
passing relatively apart from the position of the string core are apparent as faint patches outside of the string loop on the left 
panel. [Benabed & Bernardeau 2000]. 

In the case of a straight string which is aligned along the y axis, the deflection angle (or half of the deficit angle) is 
AttG^l [Vilenkin & Shcllard, 2000] and this yields a displacement £ x — ±£o with 

Co = 47rG A/ 2? lsSjS /2? s ,us (71) 

with no displacement along the y axis. I?j ss s and 2?s,us are the cosmological angular distances between the last 
scattering surface and the string, and between the string and us, respectively. They can be computed, in an Einstein- 
de Sitter universe (critical density, just dust and no A), from 

T>(z u z 2 ) = ^t^-[(1 + zi)- l/2 (1 + Z2)- 1/2 } (72) 

by taking z\ = for us and z 2 ~ 1000 for the last scattering surface; see [Bartelmann & Schneider, 2001]. For the 
usual case in which the redshift of the string z s is well below the zj sg one has 2?j gs S /T>\ ss us — l/\/l + z s . Taking 
this ratio of order 1/2 (i.e., distance from us to the last scattering surface equal to twice that from the string to the 
last scattering surface) yields z s ~ 3. Plugging in some numbers, for typical GUT strings on has Gfj, ~ 10~ 6 and 
so the typical expected displacement is about less than 10 arc seconds. Benabed & Bernardeau [2000] compute the 
resulting B component of the polarization and find that the effect is entirely due to the discontinuity induced by the 
string, being nonzero just along the string itself. This clearly limits the observability of the effect to extremely high 
resolution detectors, possibly post-Planck ones. 

The situation for circular strings is different. As shown by de Laix & Vachaspati [1996] the lens effect of such a 
string, when facing the observer, is equivalent to the one of a static linear mass distribution. Considering then a 
loop centered at the origin of the coordinate system, the displacement field can be expressed very simply: observing 
in a direction through the loop, £ has to vanish, while outside of the loop the displacement decreases as cei/a, i.e., 
inversely proportional to the angle. One then has [Benabed & Bernardeau, 2000] 



-> (X{ 

£(a) = — 2£ — a with a > aj, 



(73) 
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where ai is the loop radius. 

This ansatz for the displacement, once plugged into the above equations, yields the B field shown in both panels 
of Figure El A weak lensing effect is barely distinguishable outside the string loop, while the strong lensing of those 
photons traveling close enough to the string is the most clear signature, specially for the high resolution simulation. 
One can check that the hot and cold spots along the string profile have roughly the same size as for the polarization 
field in the absence of the string loop. The simulations performed show a clear feature in the maps, although limited 
to low resolutions this can well be confused with other secondary polarization sources. It is well known that point 
radio sources and synchrotron emission from our galaxy may contribute to the foreground [de Zotti et al. 1999] and 
are polarized at a 10 % level. Also lensing from large scale structure and dust could add to the problem. 

VI. COSMIC DEFECTS IN PERSPECTIVE 

Cosmic defects have proved very interesting and fruitful in high-energy physics and astrophysics. Their generic 
production in grand unified theories has made defects an active field of research for over two decades. Many of the 
interesting subjects now associated with defects were only briefly mentioned in these notes, like the internal structure 
of defects -leading to persistent currents in their cores- and, as a consequence, the possible generation of primordial 
magnetic fields. Also, primordial gravitational waves, extremely high-energy phenomena associated to cosmic rays 
and uhecrons, electroweak baryogenesis and, finally, the very active condensed-matter-cosmology interface, dubbed 
cosmos in the lab, equally -and unjustly- received no attention [compensate for this with references like Vilenkin & 
Shellard, 2000; Hindmarsh & Kibble, 1995; Kibble, 2002; Gangui, 2001b, for example - yes, in that order.]. With 
regards to the most transparent test of current cosmology, namely the CMB and matter power spectra, (not so) recent 
investigations have pointed out severe problems in virtually all models where cosmic defects are the main source of 
the seeds of structure in the universe 22 . In the case of cosmic strings, however, these bad news were reached by 
the use of non negligible, albeit well-founded, approximations in order to cope with the limited range of realistic 
defect simulations [bad or good news, depending on which side you are.] Although the whole method of unequal time 
correlators employed by most of the groups can be regarded as a good approximation to reality during both the matter 
and radiation eras, the important transition in between must be looked at more carefully, as the above-mentioned 
correlators do not scale as expected. Recent, full Boltzmann analyses aiming to solve this handicap are in progress 
[e.g., Landriau & Shellard, 2002, 2003] and already producing interesting results. 
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